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1 . 0  Introduction 

This  report  describes  a  stochastic  radiowave  propagation 
model  useful  for  assessing  the  effects  of  forests  upon  wide- 
bandwidth  digital  radio  communication  systems  operating  in  the 
200-2000  MHz  band.  It  represents  an  extension  of  the  model 
first  reported  in  CyberCom  Technical  Report  CTR-108-01  [42]. 

The  theoretical  basis  for  this  model  was  first  developed  by 
Foldy  [23]  and  later  extended  by  I.ax  [45],  Twersky  [68]  and  others. 
According  to  this  model,  the  forest  is  represented  as  a 
time-invariant,  stratified  configuration  of  randomly-positioned 
and  randomly-oriented  discrete  canonical  scatterers.  Tree  trunks 
are  modeled  as  inf initely-long  circular  dielectric  cylinders; 
branches  as  f initely-long  circular  dielectric  cylinders;  and 
leaves  as  flat  dielectric  discs.  The  propagating,  radiowave 
is  considered  to  consist  of  two  parts:  a  mean  (coherent) 

component  derived  by  averaging  statistically  over  an  ensemble 
of  forest  ponf igurations,  and  a  residual  random  (non-coherent , 
diffuse,  or  fluctuant)  component  uniquely  determined  by  a 
specific  forest  configuration. 

1 . 1  The  Coherent  Radiowave 

Initially,  the  emphasis  in  model  development  was  directed 
toward  characterization  of  the  mean  [coherent)  component  of 
the  propagating  radiowave.  A  physically  appealing  representation 
for  the  mean  component  was  obtained  by  postulating  that  the 
mean  field  components  satisfy  Maxwell's  equations  "in  the  mean" 
and  that  the  ensemble  of  discrete  scatterers  can  be  replaced 
by  an  equivalent  continuous  medium  described  by  an  effective 
dyadic  permittivity  e.  However,  in  contradistinction  to  earlier 
models  wherein  the  effective  permittivity  was  either  postulated 
or  deduced  from  measurement,  the  CyberCom  model  directly 
related  e  to  such  salient  biophysical  forest  parameters  as  the 
tree  trunk  number  density,  the  tree  trunk  diameter  probability 
density  function;  and  leaf  area  index. 


The  forest  model  was  refined  by  assuming  the  trees  to  be 
bounded  below  by  a  smooth  foreet  floor  and  bounded  above  by 
air.  Earlier  efforts  [15,62,641  had  suggested,  however,  that  the 
introduction  of  the  ground  complicates  the  model  significantly. 
In  the  earlier  model  [42],  these  complications  were  avoided  by 
allowing  the  forest  floor  to  recede  to  infinity  thereby  reducing 
the  stratified  forest  model  to  a  half-space  representation; 
the  newer,  more  general  model  described  in  this  report,  however, 
incorporates  a  forest  floor  of  arbitrary  permittivity.  The 
electromagnetic  boundary  value  problems  involved  in  both  models 
were  solved  using  the  classical  approach  first  described  by 
Sommerfeld  and  later  extended  by  Brekhovskikh  [5], 

The  stratified  forest  models  described  above  are 
time-harmonic  models  in  the  sense  that  the  signal  radiated  by 
the  transmitting  antenna  is  a  sinusoidal  waveform  of  angular 
frequency  ai  .  However,  because  the  equivalent  forest  continuum 
characterized  by  the  effective  dyadic  permittivity  e  is  linear, 
Fourier  transform  techniques  have  been  employed  to  generalize 
the  model  so  that  it  accommodates  arbitrarily  modulated  waveforms. 
Using  this  approach  .  •  has  been  possible  to  determine  the  mean 
pulse  distortion  and  differential  pulse  delay  associated  with 
radiowave  propagation  through  a  stratified  forest. 

1 . 2  The  Incoherent  Radiowave 

As  the  propagating  radiowave  penetrates  deeper  and  deeper 
into  the  forest,  the  relative  contribution  of  the  random 
(incoherent)  component  becomes  progressively  more  important. 

This  trend,  anticipated  from  experiments  [51,52],  was  first 
predicted  theoretically  using  a  two-dimensional,  unbounded  forest 
consisting  solely  of  tree  trunks  [42].  The  application  of  this 
early  model,  however,  is  limited  to  tree  trunk  diameters  that 
scatter  in  the  low-frequency  (Rayleigh)  regime?  for  the  200-2000 
MHz  band  the  model  proves  valid  only  for  tree  trunks  less  than 
1  centimeter  in  diameter. 

The  early  low-frequency,  Rayleigh-scattering  model  [42, 
Section  7]  has  been  improved  and  is  now  no.  longer  restricted  to 


the  model 


in  addition  to 


small-diameter  trees-  Further, 
describing  the  intensity  (power)  of  both  the  coherent  and 
non-coherent  components,  can  also  be  used  to  determine  the  space 
correlation  function  and  frequency  correlation  function  of  the 
propagating  radiowave.  Fourier  transformation  of  the  former 
describes  the  angular  spectrum  of  radiowaves  arriving  at  a  point; 
Fourier  transformation  cf  the  latter  describes  the  delay  spread 
of  the  forest  scatter  path. 

1 . 3  Scope 

This  report  describes  a  stochastic  radiowave  propagation 
model  useful  for  assessing  the  effects  of  forests  on  wide¬ 
band  digital  communication  systems.  In  Section  2  the  forest 
is  characterised,  first  as  to  the  permittivity  of  constituent 
woody  material,  then  as  to  the  geometries  of  trunks,  branches 
and  leaves.  Finally  it  is  recognized  that  forests,  in  general, 
have  canopies  supported  by  trunks,  and  some  overall  models  are 
proposed.  Section  3  presents  the  susceptibilities  and 
attenuations  for  an  unbounded  forest  of  tree  trunks  having  uniform 
or  exponentially-distributed  diameter  distributions.  Section  4 
considers  a  more  realistic  forest,  bounded  not  only  above  by 
the  air  interface  but  also  below  by  the  ground.  The  theory 
is  developed;  intra-forest  multipath  propagation  is  evaluated; 
the  concept  of  basic  transmission  loss  is  introduced;  and  Doppler 
shift  is  considered.  Section  5  presents  detailed  basic 
transmission  losses  and  pulse  responses  lor  direct,  reflected, 
dir  ect-pl,  us-reflected ,  lateral  and  total  waves  in  semi-infinite 
leaf  and  trunk  forests.  Section  6  uses  an  integro-dif f erential 
equation  for  the  space-frequency  correlation  function  to  evaluate 
the  properties  of  the  incoherent  field,  which  becomes  important 
at  the  higher  frequencies  of  interest  here.  A  two-dimensional 
forest  of  circular,  fixed-diameter,  randomly-positioned  trunks 
is  studied.  From  the  general  space-frequency  correlation  function 
are  derived  the  frequency  correlation  function  for  thin  trunks 
(yielding  coherence  bandwidth)  and  the  space  correlation  function 


for  arbitrary  diameters.  The  coherent  and  incoherent  intensities 
are  also  found  for  this  general  case,  removing  the  size 
limitations  of  the  earlier  report  [42]. 
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2 . 0  Biophysical  Forest  Characterization 

The  characterization  of  a  forest  will  be  undertaken  at 
several  levels.  It  is  first  found  in  Section  2.1  to  be  supported 
by  many  sources  that  the  basic  materials  of  green  wood  and  leaves 
(consisting  largely  of  water)  have  similar  dielectric  properties. 

A  CyberCom  model  is  established  for  current  work,  although  it 
may  be  refined  after  further  study.  The  second  level  of 

consideration  is  basically  geometrical,  but  quite  complex  in 
its  specification,  and  even  more  so  in  its  electromagnetic 

implications.  This  level  involves  the  detailed  description 
of  trunks,  branches,  and  leaves.  The  required  parameters  include 
densities,  sizes,  and  orientations,  considered  in  Section  2.2. 
The  third  level  requires  standing  back  from  the  parts  of  trees, 

and  even  from  the  trees  themselves,  to  see  the  forest,  not  as 

a  whole  it  turns  out,  but  as  a  canopy,  consisting  mostly  of 

branches  and  leaves,  supported  above  the  ground  by  the  trunks. 
This  view  is  considered  in  Section  2.3.  Finally,  forest  models 
are  synthesized  by  quantifying  the  above  concepts  in  practical 
cases.  Section  2.4  indicates  sources  of  numerical  data  and 

proposes  several  CyberCom  forest  models. 

2.1  Electrical  Properties!  of  Green  WGod  and  Leaves 

The  electrical  properties  of  green  wood  and  leaves  can 
be  specified  in  terms  of  their  relative  permittivity  and 

permeability  p^.  As  with  most  biological  materials,  the  relative 
permeability  p^  is  close  to  unity.  The  relative  permittivity 
is  complex  and,  for  the  complex  exponential  time  dependence 

exp[ju)t],  can  be  written  in  the  form 

C't  =  ”  3eJ  (2_1_1) 

where  e  ^  and  c”  represent,  respectively,  the  real  and  imaginary 
parts  of  the  permittivity  The  imaginary  part  is  proportional , 

of  course,  to  the  conductivity  o  through  the  relation 

ej  =  a/we.'  '  (2-1-2) 


v 
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where  eo  is  the  permittivity  of  free  space.  The  real  part  ££ 
is.  often  called  the  dielectric  constant;  the  imaginary  part 
e'£,  the  loss  factor. 

Although  the  electrical  properties  of  wood  have  been  the 
subject  of  numerous  studies  [10,30,37];  relatively  few  measured 
data  pertain  to  green  wood  and  leaves,  and  fewer  still  to  fre¬ 
quencies  above  100  MHz  [  38,66].  The  most  appropriate  appear 
to  be  those  of  Broadhurst  [6  ].  These  studies  suggest  that, 
firstly,  the  permittivity  of  green  wood  is  not  strongly  species- 
dependent,  although,  taken  as  a  class,  coniferous  (needle-bearing, 
tree  wood  and  needles  appear  to  have  somewhat  smaller 
permittivities  than  those  of  deciduous  (leaf-bearing)  tree  wood 
and.  leaves.  Secondly,  although  the  importance  of  intrinsic 
water  content  for  the  relative  permittivity  of  wood  has  been 
substantiated  by  many  studies,  because  green  wood  nearly  always 
has  a  high  intrinsic  water  content  (approximately  75%  by  volume), 
the  permittivity  of  green  wood  is  relatively  independent  of 
water  content.  Thirdly,  Broadhurst 's  measurements  on  Tulip  trees 
(LI,  D3  of  Fig.  2-1)  and  Bamboo  (his  Fig.  3)  suggest  very  little 
difference  between  the  permittivity  of  green  wood  and  leaves. 

Broadhurst  ‘  s  data  for  Tulip  tree  and  Maple  (L2  of  Fig.  2-1)  sug¬ 
gest  that  above  10  MHz  the  dielectric  constant  is  only  weakly 
dependent  upon  frequency  and  that  over  the  band  200-2000  MHz 
r.  ^  is  approximately  equal  to  40.  For  frequencies  below  1000 
MHz,  the  imaginary  part  e"  appears  to  be  dominated  by  conduction 

X* 

losses  and  so  decreases  linearly  with  increasing  frequency. 
Above  1000  MHz,  however,  relaxation  losses  associated  with 
molecular  polarization  begin  to  dominate  and,  in  this  frequency 
region,  begins  to  increase  with  increasing  frequency,. 

Accordingly,  in  this  study  of  UHF  radiowave  propagation  through 
forests,  the  relative  complex  permittivity  of  green  wood  and 
leaves  has  been  modeled  using  the  relation 


=  40.0 


'  j  fei 


2-0fGHz 

[1  +  {fGHz/20*0)21 


(2-1-3) 


where  f  qhz  is  the  radiowave  frequency  expressed  in  GHz. 


Table  2-1  summarizes  the  sources  used  to  determine  the  di¬ 
electric  properties  of  the  forest  constituents,  that  is,  green 
wood  and  leaves,  at  50  MHz  and  above.  Most  of  the  data  are  in 
the  form  of  dielectric  constant  and  loss,  tangent  6.  From 

these,1  the  real  and  imaginary  parts  of  susceptibility  (  X'^  nd 
X  £ ,  respectively)  have  been  found  using  the  relations  X  ^  =  ej 

-  1  and  X£  =  e  Jtan  6  .  The  results  are  plotted  in  Figures  2-1 

and  2-2.  The  CyberCom  model  is  seen  to  be  a  reasonable  approxima¬ 

tion,  but  it  may  be  noted  that,  in  general,  the  deciduous  and 
leaf  (DSL)  curves  lie  above  those  for  coniferous  and  needles 
( C&N ) .  It  seems  likely,  therefore,  that  there  may  be  a  requirement 
for  specializing  the  CyberCom  modi  1  to  deciduous  and  coniferous 
trees.  These  newer,  improved  models  may  turn  out  to  be  continuous 
functions  of  wood  density,  moisture  content,  and/or  other 
parameters. 


The  electromagnetic  properties  of  the  forest  constituents 
(trunks,  branches,  and  leaves)  depend  not  only  on  the  electrical 
properties  of  the  damp,  woody  material  considered  above,  but 
also  on  the  geometrical  factors  of  size,  shape,  and  orientation. 
The  tree  trunks  are  the  largest  forest  scatterers  and  so,  are 
considered  in  greatest  detail.  In  the  1983  CyberCom  Report  [42], 
a  forest  of  uniform  trunk  diameter  was  studied.  Here,  two  major 
forest  types  of  non-uniform  trunk  diameter  are  considered:  the 

homogeneous  forest,  characterized  by  an  even-aged  stand  with 
a  normal  (gaussian)  distribution  of  trunk  diameters,  and  the 
inhomogeneous  forest,  characterized  by  an  uneven-aged  stand  with 
an  exponential  (inverse  J-shaped)  diameter  distribution.  The 
latter  is  characterized  by  a  slope  defined  by  a  diameter  Dav 
which  may  not  be  the  true  average  diameter  because  of  tree  trunk 
thinning  by  man  or  natural  causes.  Branch  orientation  statistics 
for  conifers  have  been  found  and  are  presented.  Branch  diameter 
and  length  will  be  considered  in  a  subsequent  report.  Leaves 
are  discussed  in  some  detail  because  they  may  be  of  interest 
at  the  higher  frequencies. 
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Sources  of  Dielectric  Properties  of  Forest  Constituents 
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Figure  2-1:  Susceptibility  of  Forest  Constituents  (y 


Figure  2-2:  Susceptibility  of  Forest  Constituents  (x 


2.2.1  Trunk  Diameter  Distributions 

The  equations  presented  in  Sections  4.1  and  4.2  of  CyberCom 
Technical  Report  CTR-108-01  [42]  for  the  effective  dyadic 

susceptibility  (x)  and  specific  attenuation  (a)  of  tree  trunks 
are  predicated  upon  a  fixed  tree  trunk  diameter  and  are, 
therefore,  applicable  only  to  homogeneous  forests.  To  calculate 
a  weighted  average  of  x  and  a  for  a  non-homogeneous  forest  the 
tree  trunk  diameter  probability  density  function  (or  stand 
structure )' must  be  known. 

In  general,  the  trunk  diameter  probability  density  function 
of  an  even-aged  stand  or  plantation  is  nominally  normal  (gaussian) 
[28,54,76].  Expressed  mathematically  as 

p(D)  -  ( /2t ra)  'exp]- (D-D) 2 /2 a2}  (2-2-1) 

the  normal  probability  density  function  is  uniquely  specified 
by  two  parameters  -  the  mean  tree  trunk  diameter  (D)  and  the 

standard  deviation  (a).  The  mean  trunk  diameter  of  an  even-aged 
stand  depends  upon  the  age  of  the  stand  -  older  stands  having 
larger  diameters  ,  and  the  species  composition  -  certain  species 
grow  faster  than  others.  The  projected  average  trunk  diameters 
at  30,  60,  and  90  years  for  several  upland  deciduous  species 

are  shown  in  Table  2-2  [28].  The  standard  deviation  of  the 

trunk  diameters  about  the  mean  also  depends  upon  the  age  and 
composition  of  the  stand,  although  composition  appears  to  be 
the  more  important  parameter.  Trunk  diameter  data  taken  from 
87  half-acre  plots  used  by  the  USDA  Forest  Service  in  Ohio, 
Kentucky,  Missouri  and  Iowa  suggest  that  the  coefficient  of 
variation  (the  ratio  of  o/D)  ranges  from  about  0.7  (young, 
irregular  stands)  to  about  0.1  (old,  uniform  stands)  [refer 
to  Table  6  of  Gingrich  l 28]]. 

The  similarity  of  nearly-normal  empirical  probability  density 
functions  of  trunk  diameter  to  the  normal  probability  density 
function  can  be  measured  in  terms  of  "skewness"  and  "kurtosis." 
Skewness  is  a  measure  of  lack  of  symmetry;  a  density  function 

with  a  longer  tail  to  the  right  of  the  mean  has  positive  skewness. 
Kurtosis  is  a  measure  of  the  extent  to  which  the  height  of  the 


Table  2-2:  Projected  Average  Trunk  Diameter  (inches) 


< 

V 

Species 

30  years 

60  years 

90  years 

V 

Yellow  poplar 

6.9 

14.3 

21.9 

■ 

Black  walnut 

6.1 

12.6 

18.8 

1 

Scarlet  oak 

5.0 

11. 0 

17.9 

Red  oak 

4.6 

10.1 

16.8 

White  ash 

4.7 

10.0 

16.1 

Black  oak 

4.8 

10.1 

15.9 

v 

Sugar  maple 

3.9 

8.4 

13.5 

*2 

Beech 

3.4 

7.4 

12.1 

1 

White  oak 

3.6 

7.5 

11.8 

Hickory 

3.5 

7.1 

11.0 

v 

Chestnut  oak 

2.8 

5.8 

9.7 

Table  2-3:  Trunk  Diameter  Distributions  [Uneven-Aged  Stands] 


Author  Region  Trees  Diameter 


Prank  et  al  [24] 
Roach  [54] 

Wiant  [76] 


Northeast  US 

Allegheny 

Appalachian 


Trees 


Spruce-Fir 

hardwoods 

hardwoods 


/ 


empirical  density  function  exceeds  that  of  the  normal  density 
function;  an  empirical  density  function  with  a  sharper,  higher 
peak  than  the  normal  has  positive  kurtosis.  All  even-aged  stands 
have  some  positive  skewness,  although  skewness  decreases  with 
age  [2  ,28].  Kurtosis,  usually  positive  for  very  young  stands, 
becomes  negative  as  the  stand  ages.  Nevertheless,  both  Gingrich 
.[28]  and  Roach  [54]  have  concluded  that  the  coefficients  of 
variation,  si.ewness,  and  kurtosis  are  more  closely  related  to 
tree  average  diameter  than  to  stand  site  or  age.  This  conclusion 
could  prove  especially  helpful  in  developing  forest  models 
suitable  for  radiowave  propagation  prediction  from  forest  stand 
data.  Near ly-norma 1  probability  density  functions  for  three 
even-aged  stands  are  shown  in  Figure  2-3. 

The  Weibull  probability  density  function  has  also  been 
suggested  as  a  probabilistic  model  for  tree  trunk  diameters. 
In  [15]  an  example  is  given  for  a  20-year  old  shortleaf  pine 
(Pinus  echinata)  plantation  having  800  trees  per  acre  (1980 
trees  per  hectare).  The  approximating  Weibull  probability  density 
function,  however,  corresponds  closely  to  a  normal  probability 
density  function  having  a  mean  trunk  diameter  of  15.3  centimeters 
and  a  standard  deviation  of  3.8  centimeters. 

The  probability  density  function  of  an  uneven-aged  stand 
is  nominally  exponential  (sometimes  called  inverse  J-shaped 
in  forestry).  Expressed  mathematically  as 

p(D)  =  (1/D) exp {-D/D}  (2-2-2) 

the  exponential  probability  density  function  is  uniquely  specified 
by  a  single  parameter  -  the  mean  tree  trunk  diameter  (D).  In 
contradistinction  to  an  even-aged  forest,  the  mean  trunk  diameter 
of  a  mature  uneven-aged  stand  does  not  depend  upon  its  age  - 
only  upon  its  species  composition.  Tabie  2-3  shows  several 
sources  of  uneven-aged  forest  data  (  one  of  which  is  plotted 
in  Figure  2-4  )  from  the  northeastern  United  States,  and  the 
mean  trunk  diameters  derived  by  fitting  exponential  probability 
density  functions  to  the  data  using  a  least-squares  method. 
The  mean  trunk  diameters  are  all  close  to  2.5  inches  (6.35  cm); 
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Figure  2-4:  Trunk  Diameter  Probability  Density  Function 

I  Uneven-aged  Stand]  [76] 
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the  corresponding  correlation  coefficients  of  the  least-squares 
fit  are  all  near  unity  and  suggest  a  very  good  fit. 

It  should  be  realized,  however,  that  the  parameter,  mean 
tree  trunk  diameter  (D),  as  used  above  in  Equation  (2-2-2)  is 
strictly  a  mathematical  measure  for  the  slope  of  the  exponential 
probability  density  function  and  that  the  mean  tree  trunk  diameter 
as  measured  in  the  field  may  be  considerably  larger,  especially 
if  the  smaller  trees  are  not  well  represented  because  of  thinning 
and/or  competition  by  the  larger  trees. 

The  previous  paragraphs  underscore  the  importance  of  mean 
trunk  diameter  in  characterizing  the  trunk  diameter  probability 
density  function.  Other  factors  bearing  on  the  density  function, 
while  not  so  important  as  mean  trunk  diameter,  are  the  type 
of  stand  (even-aged  or  uneven-aged)  and,  if  even-aged,  its  age 
and  species  composition  (especially  the  latter).  In  addition 
to  the  trunk  diameter  probability  density  function,  the  trunk 
number  density  is  also  required  to  estimate  the  mean  effective 
dyadic  susceptibility  and  specific  attenuation  of  non-homogeneous 
forests.  A  potentially  useful  aid  in  characterizing  forests 
for  radiowave  propagation  prediction  may  be  the  forester's 
"stocking  guide."  These  guides  relate  trunk,  (stem)  number 
density,  trunk  basal  area,  mean  trunk  diameter  in  the  manner 
exemplified  by  Figure  2-5.  Silviculturists  use  these  guides 
to  manage  forest  yield.  As  noted  by  Roach  ,[54],  "the  guide 
has  a  rational  biological  basis;  it  is  generally  independent 
of  the  influence  of  site,  age,  or  stand  structure;  and  the 
measurements  that  are  needed  for  its  application  can  be  easily 
and  quickly  obtained  in  the  field."  The  trunk  number  density 
and  basal  area  have  been  summarized  by  Cannell  [  9  ]  for  many 
forests  of  the  world;  the  corresponding  mean  trunk  diameter 
and,  perhaps,  species  composition  can  be  inferred  from  stocking 
guides;  and  from  these  data  the  trunk  probability  density  function 
can  be  estimated. 

2.2.2  Branch  Orientation  , 

The  equations  presented  in  Sections  4.1  and  4.2  of  Cyber Com 
Technical  Report  CTR-108-01  for  the  effective  dyadic 
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susceptibility  and  specific  attenuation  of  branches  and  leaves 
involve  the  probability  density  functions  of  their  inclination 
angles.  Probability  density  functions  for  the  branch  inclination 
angles  of  80-year  old  Lodgepole  pine  in  Wyoming  were  measured 
by  Gary  [27]  and  are  presented  here  in  Figure  2-6(a).  It  is 

apparent  from  this  figure  that  the  inclination  angles  of  the 

younger  branches  in  the  upper  crown  are  steeper  than  those  of 
the  older  branches  in  the  lower  crown.  This  trend  is  confirmed 

by  Ford  [24]  and  by  Kimes  et  al  [39].  The  negative  angles 

indicate  that  some  branches  pointed  downward.  Gary  [27]  also 
reported  that  there  was  no  significant  departure  from  azimuthal 
symmetry.  Linear  fits  to  these  data  have  been  made  by  CyberCom 
and  are  presented  in  Figure  2-6 (b).  These  linear  models  for 
the  probability  density  function  of  the  branch  inclination  angle 
will  be  employed  by  CyberCom  in  subsequent  numerical  evaluations 
of  the  effective  dyadic  susceptibility  and  specific  attenuation 
of  coniferous  forests.  Analogous  data  for  the  branch  inclination 
angles  of  deciduous  trees  have  not  yet  been  located. 

2.2.3  Leaf  Properties 

Because  leaf  sizes  approach  the  wavelength  only  at  the 
higher  frequencies  of  interest  (X  at  2000  MHz  =  6"),  it  appears 
that  the  only  leaf  parameter  of  real  importance  in  the  present 
model  is  the  fractional  volume  occupied.  This  is  considered 
in  the  overall  forest  characterization  of  Section  2.4.  However 
because  of  the  previous  concentration  on  leaf  parameters  in 
CyberCom  Technical  Report  CTR-108-01  and  because  of  possible 
future  interest  (at  higher  frequencies  for  the  applications), 
the  results  of  a  study  of  leaf  parameters  is  included  below. 

2. 2. 3.1  Leaf  Diameter  and  Thickness 

Two  good  sources  for  the  area  and  thickness  of  leaves 
of  southeastern  deciduous  trees  were  found.  The  data  from 
Rothacher  [56]  are  summarized  in  Table  2-4 The  composition 
by  species  for  the  Tennessee  Valley  shows,  in  the  weight  column 
and  in  the  numbers  of  trees  larger  than  11  inches  diameter  at 
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Figure  2-6 (a):  Branch  Inclination  Angles  (measured) 
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Figure  2-6 (b):  Branch  Inclination  Angles  (CyberCom) 
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breast  height  (dbh),  the  predominance  of  oaks.  The  numerous 
small  sourwoods  form  an  understory  following  a  heavy  cut  15 
years  before  data  acquisition.  Dividing  the  total  leaf  areas 
by  the  numbers  of  leaves  according  to  species,  the  area  per 
leaf  column  is  obtained  with  areas  ranging  from  2.5  cm2  for 
pine  needles  to  159  cm2  for  hickory  leaves.  Catalpa  leaf  areas 
up  to  7 _>0  cm2  have  been  reported  11),  but  the  more  representative 
values  of  about  55  for  Scarlet  Oak  and  Yellow  Poplar  give  a 
typical  equivalent  circular  diameter  of  about  8  cm.  The  average 
leaf  thicknesses,  from  Carpenter  and  Smith  [111,  are  for  shade 
leaves  in  Lexington,  Kentucky.  They  range  from  0.15  mm  to  0.30 
mm,  with  a  typical  value  of  0.2  mm.  The  volumes  per  leaf  range 
from  0.2  to  2.9  cc.  The  specific  gravities  are  derived  by 
dividing  the  total  leaf  weights  in  grams  by  the  total  volumes 
in  cubic  centimeters.  The  values  range  from  0.14  g/cc  to  0.57 
g/cc  with  a  typical  figure  of  about  0.35  g/cc. 

The  data  for  conifers  provided  in  Table  2-5  are  from  the 
Biological  Handbook  [  1  ]  except  the  first  two  items  from  Rothacher 
[56].  Needle  lengths  have  been  added  from  Little  [46].  An 
effective  cylinder  diameter  for  the  short-leaf  pine  needle  may 
be  obtained  as  0.5  mm  using  a  length  of  20  cm  to  account  for 
two-leaf  fascicles  per  Rothacher. 

2. 2. 3. 2  Foliage  Inclination  Angles 

Probability  density  functions  for  the  needle  inclination, 
angles  of  Lodgepole  pine  were  measured  by  Kimes  et  al  [39]  and 
found  to  be  relatively  invariant  with  height.  CyberCom  has 
found  that  these  density  functions  can  be  approximated  by 

!l/60 ,  (O°<0<3O#) 

(2-2-3-2-1) 

(90-0) /3600 ,  (30°<0<90*) 

Subsequent  measurements  on  Douglas  fir  by  Smith  et  al  [60] 
substantiate  these  results.  Smith  et  al  also  considered  the 
leaf  inclination  angles  in  a  deciduous  oak-hickory  canopy.  Here, 
only  the  lower-third  of  the  canopy  agreed  with  the  coniferous 
probability  density;  the  two  upper-thirds  had  probability  density 
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Table  2-5  Total  Surface  Areas  of  Conifer  Needles  [  1) 


» 


r 


functions  that  wore  essentially  uniform  for  inclination  angles 
between  0  and  33  degrees  from  the  horizontal. 


Forest  Canor 


Structure 


The  most  obvious  structural  feature  of  a  forest  is  not 
just  that  it  Consists  of  trunks,  branches  and  leaves,  but  that 
these  are  formed  into  a  canopy  of  branches  and  leaves,  supported 
by  the  trunks.  (For  isolated  or  widely-spaced  trees,  the 
continuous  canopy  reduces  to  an  ensemble  of  crowns) ;  This  simple 
concept,  with  importance  for  electromagnetic  modelling,  is  not 
well  covered  in  the  literature.  The  lumberman  pays  for  wood 
without  regard  to  its  source  and  the  business  expert  cares  no 
more  about  the  growing  heights  of  leaves  or  branches.  Only 
academic  sources,  sometimes  in  support  of  the  above 
economically-oriented  workers,  were  found  for  canopy  structure, 
largely  coniferous.  Some  items  considered  below  include  foliage 
distribution  and  canopy  thickness. 


The  distribution  of  foliage  weight  with  height  in  the  canopy 
of  an  even-aged  stand  (e.g.,  a  plantation  of  trees  planted  at 

the  same  time)  can  be  reasonably  well  described  by  a  normal 
(gaussian)  probability  density  function  [24,27,63,73]. 
Nevertheless,  as  might  be  expected,  not  all  even-aged  stands 
show  good  agreement.  For  example,  although  Waring  et  al  [67] 
showed  this  model  to  be  valid  for  Ponderosa  pine,  they  found 

Douglas  fir  and  Grand  fir  had  modal  peaks  well  below  mid-canopy 
height.  Further,  although  Stephens  [63]  found  this  model  suitable 
for  Red  pine  and  Japanese  cedar,  he  also  found  the  weight 

distributions  for  three  species  of  deciduous  trees  from  Japan 

to  have  modal  peaks  well  above  mid-canopy.  This  upward  skew 
of  the  weight  distribution  for  deciduous  trees  is  also  supported 
by  Figure  65  of  Assmann  (  2  ]  and  by  Table  2-6  from  Smith  et 
al  (601. 


Although  canopy  thickness  data  are  scarce,  Ford’s  curve 
(241  for  a  stand  of  Sit.a  spruce  suggests  a  ratio  of  canopy  thick¬ 
ness  to  forest  height  of  about  0.6.  Gary  (27],  too,  suggests 
a  ratio  close  to  0.6  for  Lodgepole  pine;  Kinerson  and  Fritschen 
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Table  2-6:  Canopy  Layer  Heights  and  Leaf  Area  Index  (LAI) 


Canopy 

Layer 

Douglas 
Median 
Height  (m) 

Fir 

LAI 

Oak-Hickory 

Median 

Height  (m)  LAI 

Upper 

23.2 

1.5 

18.3 

3.4 

Middle 

14.0 

5.3 

11.0 

CO 

• 

o 

Lower 

4.7 

1.0 

3.7 

0.4 

Note:  LAI  -  Total  one-sided  leaf  area  /  projected  area  under  canopy 


Table  2-7:  Total  Green  Weight  of  Above 


-Ground  Tree  Biomass 


(millions 

of  metric 

tons)  ( 

excluding  foliage)  1 

[48] 

Softwoods  (Coniferous) 

Ha 

rdwoods  (Deciduous) 

Southern  Pines 

2175 

35.0% 

Oa 

ks 

1543 

25% 

Loblolly 

923 

15.0% 

So 

ft  Hardwoods 

1542 

25% 

Slash 

461 

7.5% 

Other  Pines 

791 

12.5% 

•  , 

Other  Softwoods 

281 

4.0% 

Ot 

her  Hardwoods 

717 

11% 

.0% 


Total  Softwoods  2455  39 


Total  Hardwoods  3803  61% 


[40]  about  0.5  for  Douglas  fir;  and  Assmann  [  2  }  about  0.4  for 
Norway  spruce  and  0.3  for  Scots  pine.  All  are,  obviously, 
conifers.  For  deciduous  trees,  Hutchinson  and  Matt  [33]  give 
a  relative  canopy  thickness  of  0.5  for  Yellow  Poplar  at  Oak 
Ridge,  Tennessee,  and  Assmann  [  2  ]  gives  a  value  of  0.5  for 
Common  beech.  In  view  of  these  data.  Cyber Com  will  employ  a 
value  of  0.5  for  the  ratio  of  canopy  thickness  to  forest  height 
for  all  forest  models,  coniferous  and  deciduous. 

2.4  Forest  Models 

Although  models  of  all  forests  of  the  world  may  ultimately 
be  of  interest,  it  seems  best  to  start  with  a  region  of  high 

interest  and  high  data  availability,  such  as  the  southeastern 
United  States.  Recently  empirical  test  data  were  acquired  in 
Florida  and  Tennessee  in  areas  dominated  by  oaks  and  pines. 
Table  2-7,  which  refers  to  commercial  forest  land  in  the  seaboard 
states  from  Virginia  to  Florida  shows  that  pines  comprise  35% 
and  oaks  25%  of  the  above  ground  forest  biomass. 

The  best  forest  data  base  available  -  broadest  in  both 

geographical  coverage  and  in  parameters  considered  -  appears 
to  bo  a  1982  compilation  of  "World  Forest  Biomass  and  Primary 

Production  Data"  by  M.G.R.  Cannell  [9  ].  The  first  half  of 
Table  2-8  summarizes  12  pages  of  this  book,  including  most  of 
the  relatively  complete  data  from  the  southeastern  United  States, 
largely  from  the  Great  Smoky  Mountains.  From  these  data  are 
derived  the  fractional  volumes  in  the  right  hand  columns  and 
several  representative  parameters.' 

The  ,  first  numerical  column  of  Table  2-8  gives  the  number 
of  trees  generally  larger  than  2  cm  in  diameter  at  breast  height 
(dbh)  per  hectare  (104m2)  of  forest.  A  typical  figure  of  2500 
corresponds  to  an  average  trunk  density  Pt  of  0.25  trunks/m2 
or  an  average  spacing  between  trunks  of  about  2  meters.  The 
average  tree  heights  presented  in  the  second  column  are  used 
in  establishing  the  height  of  the  forest  and  estimating  leaf 

and  branch  fractional  volumes.  In  the  third  column,  basal  area 
gives  the  total  of  the  cross-sectional  areas  of  all  the  trees 
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Table  2-8  Southeastern  United  States  Forest  Parameter 


U-f 

05 

X 

O 

<D 

T! 

C 

•4 

< 

M 

C 

O 

•H 

4J  ca 
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U  -U  <TJ 
[14  to  P 
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4J 

co 
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C 
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•r^ 

CO 

w 

w 

.c 

s 

u 

« 

e  jc 

<U  V4 

>4 

■U  (13 

Q 

t/5  CQ 

(dj^2cm)  per  hectare.  A  typical  basal  area  of  45  m2/ha  with 
1  500  trees/ha  gives  a  cross-sectional  area  per  tree  of  0.03  m2, 
corresponding  to  an  average  tree  diameter  of  about  20  cm.  The 
fourth  column,  stem  volume,  is  the  total'  of  all  trunk  volumes 
independent  of  :lbh. 

Dry  masses  are  shown  next  in  columns  five,  six  and  seven 
for  stem  wood  and  bark  combined,  branches,  and  leaves.  An  arrow 
indicates  that  a  figure  includes  that  for  the  adjacent  column. 
It  is  of  interest  to  note  that  the  maximum/minimum  ratio  for 
trunks  and  branches  runs  near  11,  but  for  leaves  only  1.5.  This 
may  reflect  a  basic  requirement  for  efficient  capture  of  radiation 
by  the  ennopv.  Specific  gravities  (green  vol./dry  wt.)  from 
the  Wood  Handbook  (70]  are  listed  next  in  column  eight.  These 
are  used  to  estimate  fractional  volumes  for  use  in  equations 
for  specific  attenuation.  Several  of  these  values,  notably 
those  for  maple  and  the  oaks,  differ  significantly  from  those 
calculated  in  Table  2-4  for  dry  leaves.  The  difference  may 
arise  from  the  use  of  green  leaf  thicknesses  in  the  calculations, 
but  it  is  not  understood  in  that  case  why  the  values  for  Yellow 
Poplar  arc  not  equally  far  off. 

The  weights  per  hectare  divided  by  the  specific  gravities 
give  the  volume  per  hectare,  but  to  convert  this  to  fractional 
volume  requires  the  "canopy  volume",  that  is,  its  effective 
thickness  over  the  hectare.  Although  illustrations  in  tree 
manuals  suggest  that  the  leaves  are  confined  to  about  eight-tenths 
of  the  average  tree  height,  pictures  taken  within  forests  reveal 
that  leaves  or  needles  may  be  confined  to  much  smaller  fractional 
height  and  so  one-half  was  used.  These  results  are  expressed 
in  percent  for  compactness,  and  show  high  variability. 

The  leaf  area  index  is  the  total  one-sided  leaf  area  divided 
by  the  projected  area  under  the  leaves.  Thus,  the  total  number 

of  leaves  per  hectare  is  the  product  of  the  leaf  area  index 

and  the  projected  area,  104m2,  divided  by  the  one-sided  area 
of  a  sinqle  leaf.  The  loaf  number  iensity  can  be  obtained 

by  dividing  the  number  of  leaves  per.  hectare  by  the  "canopy 

volume",  which  is  projected  area  times  effective  canopy  thickness. 
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For  example,  for  a  forest  with  a  leaf  area  index  of  6,  an 
effective  canopy  thickness  of  16  meters  and  a  per-leaf  area 
of  50  cm2,  p£  is  about  75  leaves  per  cubic  meter. 


3.0  Unbounded  Forest  Susceptibility  and  Specific  Attenuation 


This  section  of  this  report  will  investigate  the  effects 
of  changing  forest  parameters  on  the  effective  susceptibility 
and  specific  attenuation  of  the  forest  considered  as  a  uniform 
dielectric  medium.  Three  main  types  of  forest  are  considered. 
Homogeneous  (even-aged)  forests  of  tree  trunks  of  different 
diameters  from  1  to  16  inches  diameter  are  covered  in  Section 
3.1.  Non-homogeneous  (uneven-aged)  forests  of  trunks  with 
diameters  following  the  exponential  distribution  are  covered 
in  Section  3.2.  Finally,  forests  of  leaves  are  covered  in  Section 
3.3. 

3 . 1  Homogeneous  Forest  of  Tree  Trunks 

The  sensitivity  to  tree  trunk  diameter  of  the  effective 
dyadic  susceptibility  (x)  and  the  specific  attenuation  (a)  of 
an  unbounded,  homogeneous  forest  of  parallel  trees  is  shown 
in  Figures  3A-1  to  3A-18.  These  results,  which  encompass  a 
range  of  trunk  radii  extending  from  0.0125  meter  to  0.20  meter, 
have  been  calculated  using  Equation  (4-1-1-1)  and  (4-2-1-3) 
of  CyberCom  Technical  Report  CTR-108-01.  Because  of  the 
relatively  weak  dependence  of  on  frequency,  both  x  and  a  are 

sensitive  to  frequency  and  trunk  radius  primarily  through  the 
parameter  k0a.  As  a  consequence,  the  curves  presented  in  these 
figures  are  very  nearly  frequency-shifted  replicas  of  each  other. 


rt 

although  they 

are 

displaced 

in  magnitude  by 

an  amount  proportional 

to  a  power 

of 

a .  In 

the  case  of  x, 

the 

displacement  is 

proportional 

to 

a2;  in 

the  case  of  a, 

the 

displacement  is 

proportional 

to  . 

a2  in  the 

low-frequency 

(Rayleigh)  region  and 

c 

to  a  J/4  in 

the 

resonant  region.  At 

the 

higher  resonant 

frequencies. 

the 

resonant 

maxima  of  (x£» 

Xz) 

decrease  roughly 

as  the  square 

of 

the  frequency;  conversely,  the 

resonant  minima 

of  (x£,x£)  decrease  roughly  linearly  with  increasing  frequency. 
The  first  null  in  Xz  is  due  to  a  change  in  sign  of  Xz  occurring 
as  a  consequence  of  resonance  when  the  circumference  of  the 
tree  equals  the  wavelength  within  it  (300/fMHZ/e).  This  resonant 
behavior  may  be  attributed  to  circumferential  surface  waves 
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that  reinforce  themselves  after  successive  trips  around  the 
perimeter  [49,50].  It  is  apparent  from  Figures  3B-3  through 
3B-5  that  the  specific  attenuations  for  both  horizontally  and 
vertically  polarized  radiowaves  have  nearly  the  same 
high-frequency  asymptote. 

3 . 2  Non-homogeneous  Forest  of  Tree  Trunks 

In  order  to  determine  the  effective  dyadic  susceptibility 
and  specific  attenuation  in  a  non-homogeneous  forest,  wherein 
the  tree  trunk  radius  a  is  a  random  variable  described  only 
statistically  by  some  probability  density  function  p(a),  it 
is  necessary  to  effect  a  weighted  average  of  the  X  and  a  attri¬ 
butable  to  some  particular  radius  a.  In  a  non-homogeneous  forest 
the  mean  effective  dyadic  susceptibility  and  the  mean  specific 
attenuation  a  can  be  determined  from  the  following  relations: 


X  =  J x  (a)  p  (a)  da 
ci  =  J* a  (a)  p  (a)  da 


(3-1-1) 

(3-1-2) 


where , 


X(a)  =  effective  dyadic  susceptibility  for  trunk  of  radius  a 
a  (a!)  =  specific  attenuation  for  trunk  of  radius  a 


p(a 

The 

meters ) 


)  =  probability  density  function  of  trunk  radii. 


probability  density  function  of  the  trunk  radii  (dia- 
has  been  discussed  and  described  in  Subsection  2.2.1 
and  fourid  to  be  nearly  normal  (gaussian)  for  even-aged  forests 
and  exponential  for  uneven-aged  forests.  However,  in  th 
remainder  of  this  subsection,  only  uneven-aged,  rton-homogeneous 
forests  will  be  considered.  In  consonance  with  Section  2.2.1 
and  Tablb  2-4  it  will  be  assumed  that 


where  a 

The 

attenuat 


p(a)  =  ( 1/a )exp[ -a/a ] 

f  3.2  cm  is  the  mean  tree  trunk  radius. 


(3-1-3) 


mean  effective  dyadic  susceptibility  and  mean  specific 
lion  in  an  unbounded,  uneven-aged,  non-homogeneous,  trunk- 
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dominated  forest,  wherein  tree  trunk  radii  are  described  statis¬ 
tically  by  the  exponential  probability  density  function  of 
Equation  (3-1-3),  has  been  evaluated  using  the  approximations 

_  N 

X  ~  ^JCfa^p  (ai)  Aa.^  (3-1-4) 

i 

N 

a  -  (ai)p  (ai)  Aai  (3-1-5) 

i 

and  plotted  in  Figures  3-1  and  3-2.  In  effecting  these  calcula¬ 
tions  N  was  set  equal  to  12  and  the  a^  assumed  the  values  from 
1  to  12  inches.  For  convenience,  these  results  have  been 
normalized  to  a  tree  trunk  number  density  of  p  =  o.l  tree  trunk 
per  square  meter.  The  usual  smoothing  effect  of  integration 
is  very  apparent. 

In  any  real  forest,  the  exponential  probability  density 
function  proves  valid  only  over  a  range  of  a  bounded  by,  say, 
amin  to  amax*  Outside  of  this  range,  the  exponential  probability 
density  function  usually  over-estimates  the  relative  number 

of  tree  trunks.  Recognizing  this,  the  validity  of  Equations 

(3-1-4)  and  (3-1-5)  if  the  limits  of  integration  are  restricted 
to  the  range  [amj.n,amax]  can  be  assessed  by  considering  the 
normalized  distribution  function 

/amax  ~°° 

f (a) p (a) da/ j  f(a)p(a)da  (3-1-6) 

0  0 

where  f(a),  assumes  the  role  of  either  X  or  a-  More  specifically, 
if  fta)  =  X  or  if  f(a)  =  a  and  attention  is  restricted  to  the 
low-frequency  (Rayleigh)  region,  then  f(a)  «  a2  and 

flW  *  11  '  (1  +  amax/?+aLx/2j2)axP{-ar,ax/?)1  ,3*l'7) 
If  f(a)  =  a  and  attention  is  restricted  to  the  resonant  region, 
then  f(a)  *  a2/4  and 

F(amax )  *  Y ( 1 . 75 ,amax/a )/T (1 • 75)  (3—1—8) 

where  Y(x,y)  is  the  incomplete  gamma-function  defined  by 


Susceptibility  *  X 
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Figure  3-2:  Specific  Attenuation  (Trunks,  Uneven-aged 


(3-1-9) 


e 


Y(x,y)  =  /  tx  ^xpt-tldt 
0 

and  r ( x )  is  the  complete  gamma-function  defined  by 

/CO 

tx_1exp{-t}dt  (3-1-10) 

0 

The  normalized  distribution  function  F(amax)  has  been  evaluated 
in  both  the  Rayleigh  and  resonant  ragions  using  Equations  (3-1-7) 
and  (3-1-8),  respectively,  and  plotted  in  Figure  3-3.  It  is 
apparent  from  this  figure  that  in  an  uneven-aged  forest  typical 
of  the  northeastern  United  States  about  80%  of  the  mean  effective 
dyadic  susceptibility  or  the  mean  specific  attenuation  in  the 
Rayleigh  region  can  be  attributed  to  trees  having  trunk  radii 
in  the  range  4-16  centimeters  (diameter  at  breast  height  in 
the  range  3-13  inches).  About  80%  of  the  specific  attenuation 
in  the  resonant  region  can  be  attributed  to  trees  having  trunk 
radii  in  the  range  1-11  centimeters  (diameter  at  breast  height 
in  the  range  1-9  inches).  In  the  frequency  range  200-2000  MHz 
where  all  trees  having  trunk  radii  greater  than  about  1  centimeter 
(2  inch  diameter  at  breast  height)  are  resonant,  the  mean  specific 
attenuation  is  primarily  attributable  to  tree  trunk  radii  in 
the  range  1-11  centimeters  (1-9  inches  at  breast  height ) . 

3 . 3  Forest  of  Leaves 

Figures  3-4  and  3-5  show  similar  x  and  a  plots  for  the 
leaf  forest  considered  in  CyberCom  Technical  Report  CTR-108-01 
in  order  to  make  this  report  complete  in  itself  and  to  present 
comparisons  between  leaf  and  trunk  forest  characteristics*. 
Because  the  leaves  are  small  compared  to  the  wavelength  (Rayleigh 
or  hen-resonant )  the  plots  are  relatively  smooth.  The  plots 
for  large  diameter  trees  are  fairly  smooth  due  to  averaging, 
resulting  from  the  large  number  of  wavelengths  on  a  circumference. 
Because  the  trunks  are  all  vertical  and  the  leaves  are  largely 
horizontal  (0°  to  30®  tilt),  X  t>)(t  for  leaves  and  vice  versa 

•Figures  4-1,  2,  4,  6,  and  9  of  Cyber  Cow  Technical  Report  CTR-108-01  arb 

shown  as  Figures  3A-4,  3-4,  3B-4,  and  3-5. 
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Figure  3-5:  Specific  Attenuation. {Leaves,  FV  *  .1%) 
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for  trunks.  It  may  be  seen,  of  course,  that  the  a  curves  are 
closely  related  to  those  for  X"  (imaginary  plot  of 
susceptibility). 

The  sensitivity  of  X"  for  a  non-homogeneous  forest  (1-12-1) 
to  addition  of  leaves  was  briefly  investigated.  For  the  most 
interesting  case  with  vertical  antennas,  x"z  showed  only  a  small 
change  at  highest  frequencies  upon  addition  of  0.1%  fractional 
volume  of  leaves,  as  in  CyberCom  Technical  Report  CTR-108-01. 

I» 

In  the  other  case,  Xt  showed  small  changes  at  low  frequencies 
with  a  leaf  fractional  volume  of  0.003%. 
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ANNEX  3B 

Specific  Attenuation 

^  of  Homogeneous  Trunk  Forests 


Radius 
(meters ) 

Diameter 
(inches ) 

Figure 

Number 
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Figure  3A-5(b):  Effective  Dyadic  Susceptibility 
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4 . 0  Anisotropic  Forest  Slab  Model 

In  this  section  the  medium  will  be  bounded  to  represent 
the  forest.  In  l 42]  the  upper  air-forest  interface  was  introduced 
to  produce  the  half-space  model,  shown  in  Figure  4-1 (a).  Although 
relatively  simple,  it  gives  rise  to  the  lateral-wave  which  is 
essential  to  long-range  transmission.  The  slab  model  to  be 
considered  here  adds  the  lower  ground  interface  as  in  Figure 
4-1 (b).  Further  possible  model  refinements  are  the  multi-slab 
model,  representing  the  canopy  as  in  Figure  4-l(c),  and  the 
inclusion  of  height-dependent  properties  as  in  Figure  4-1 (d). 

Section  4.1  formulates  the  problem,  introducing  the 
reflection  coefficient  at  the  ground  as  well  as  at  the  air 
interface.  Section  4.2  performs  the  asymptotic  (approximate) 
evaluation  of  the  equation  just  developed  using  methods  of  [42], 
but  these  are  complicated  by  multiple  reflection  of  both  space 
and  lateral  waves  at  the  upper  and  lower  forest  boundaries. 
This  multipath  nature  is  treated  in  detail  in  Section  4.3.  Curves 
of  reflection  coefficient  versus  distance  for  the  media  of 
interest  are  shown.  Interference  of  direct  and  reflected  waves 
is  expanded  to  include  the  many  possible  cases  of  multiple 
reflection  above  and  below  the  forest.  Then  the  lateral  wave 
is  included.  This  may  be  excited  after,  or  received  after, 
multiple  reflections.  Section  4.4  introduces  the  concept  of 
basic  transmission  loss,  the  normal  tool  of  radio  link  design. 
It  is  useful  in  the  broad-band  application  to  spread-spectrum 
transmission  of  interest  here.  Finally  Section  4-5  finds  that 
terminal  motion  produces  negligible  Doppler  shifts. 

4.1  Model  Formulation 

The  anisotropic  slab  model  for  UHF  radiowave  propagation 
through  forests  is  shown  in  Figure  4-2.  The  electrically  ani¬ 
sotropic  forest  slab,  sandwiched  between  the  electrically  iso¬ 
tropic  air  and  ground,  is  characterized  by  the  effective  dyadic 
permittivity  ee0  (at  least  so  far  as  the  mean  fields  are  con¬ 
cerned):  The  scalar  permittivities  of  the  air  and  the  ground 
are  denoted  by  e0  and  £ge0,  respectively.  All  three  media  are 
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Figure  4-1:  Stratified  Forest  Models 
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Figure  4-2:  Slab  Model  for  the  Forest 


assumed  to  be  magnetically  isotropic  and  characterized  by  the 
free-space  permeability  uo .  The  transmitting  antenna  is  repre¬ 
sented  as  a  vertical  electric  (Hertzian)  dipole  having  a  time- 
harmonic  current  moment  Idi-exp  { ja> t }  located  at  height  zo  above 
the  ground. 

In  an  unbounded,  anisotropic  forest  characterized  by  a 
uniaxial  effective  dyadic  permittivity 


*  +  +  ez££ 


(4-1-1) 


the  mean  electromagnetic  fields  associated  with  a  z-directed 
(vertical)  Hertzian  dipole  can  be  derived  from  a  Hertz  vector 
potential  A  that  has  only  a  z-directed  component,  viz. 


Az(r,z)  = 


-IdS.  f  exP{"j  (J?t*£+T  I  z~z»  I )  ) 


(4-1-2) 


8ir2u)e0e. 


where 


t  =  x »  -  j  |  t  "  1  =  /et/e2  AS  ez-62 


(4-1-3) 


The  electromagnetic  fields  can  be  derived  from  the  vector 
potential  A  using  the  relations 


E  *  V (V*A)  +  kSe^A 


(4-1-4) 


H  =  ju>e,et(?xA) 


(4-1-5) 


For  the  anisotropic  forest  slab  model,  the  tangential  com¬ 
ponents  of  the  electromagnetic  fields  must  be  continuous  across 
the  air-forest  and  the  forest-ground  interfaces.  If  the  source 
of  these  fields  is  a  z-directed  (vertical).  Hertzian  dipole, 
the  boundary  conditions  can  be  succinctly  summarized  in  terms 
of  a  Hertz  vector  potential  having  only  a  z-directed  component 
by  writing 


Az(r,H)  =  et  A*J(r,H) 


etAzI(-,°)  =  egAzII(£'0) 


(4-1-6) 


(4-1-7) 
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3 [A*(r,H) ]/3z  =  3[A^I(r,H) ]/3z 

2  Z 


(4-1-8) 


3[A“(r,0)l/3z  =  3[A*XI(r,0)  ]/3z 


(4-1-9) 


where  the  superscripts  I,  II,  III  identify,  respectively,  the 
Hertz  potentials  intrinsic  to  the  air,  forest,  and  ground. 

As  a  consequence  of  these  boundary  conditions,  the  radiation 
condition  and  Equation  (4-1-2),  CyberCom  has  found  the  Hertz 
potential  within  the  anisotropic  forest  slab  to  be* 


*11 ,  X 

Az  (r,z)  = 


-Id& 
8rr2uje  .e 


[l+Rgexp{-j2T2ze}]  [l+Raexp{-j2T2 (H-z) }] 


[ l-RaRgexp{ - j 2t 2H } ] 


exp{-j (3t*r+x2 | z— z „ | ) } 


(4-1-10) 


where, 


T2  ~  GtTl 

Ra<£t>  ■ 

C  ,T«  -  C.T, 

Ra(6t)  =  r 2~r-F~ 
g  c  egT2  etx3 


(4-1-11) 


(4-1-12) 


are,  respectively,  the  Fresnel  reflection  coefficients  associated 
with  the  air-forest  and  forest- ground  interfaces,  and  where 


Ti  =  T1  ”  ^  I T £  I  3  ^  -  ej 


(4-1-13) 


x2  =  T2  "  j  I  T£  I 


'et/eg  /k«ez"6t 


(4-1-14) 


X3  =  T*  -  j  |  T  ”  |  =  AUq  - 


(4-1-15) 


are  the  dispersion  relations.  The  electromagnetic  fields  within 
the  anisotropic  forest  slab  can  be  derived  by  substituting 
Equation  (4-1-10)  into  Equations  (4-1-4)  and  (4-1-5). 

Equations  (4-1-5),  (4-1-5)  and  (4-1-10)  constitute  the 
formal  solution  for  the  electromagnetic  fields  of  a  vertical 
electric  (Hertzian)  dipole  in  an  anisotropic  forest  slab. 


‘strictly  speaking,  Equation  (4-1-10)  is  valid  oriiy  when  z>z*.  If  z<z»  , 
then  z  and  z  must  be  interchanged. 
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Asymptotic  Evaluation 

The  integral  representation  afforded  by  Equation  (4-1-10) 
for  the  Hertz  potential  within  the  anisotropic  forest  slab  is 
not  amenable  to  exact  analytic  evaluation.  However,  analytic 
asymptotic  approximations  have  been  derived  based  upon  CyberCom 1 s 
earlier  development  of  an  anisotropic  forest  half-space  model 
f 4 2  1  and  the  pioneering  efforts  of  Sachs  and  Wyatt  [57.],,  Tamir 
[65,64,19],  and  others  on  the  isotropic  forest  slab  model. 

Before  proceeding  with  the  development  of  these  asymptotic 
approximations,  it  will  prove  expedient  to  recast  the  expression 
for  the  Hertz  potential  within  an  anisotropic  forest  half-space 
[42,  Equation  (5-2-39)]  in  geometric  terms  consistent  with  Fig. 
4-1.  Thus,  introducing  the  following  correspondences  apparent 
by  comparison  of  Fig.  4-2  and  Fig.  5-2  of  [42], 

d-*(H  -  Zo)  (4-2-1) 


z  { H  -  z ) 


(4-2-2) 


into  the  equation  previously  cited,  the  Hertz  potential  within 


the  anisotropic  forest  slab  AjJ 


may  be  written  in  the  form* 


Az  (r , z) 


-Id& 
8n2d je0e. 


— 

/exp{-jx2|z  -z.| } 
r2 


+  W 


exp{-jT2 (2H-z-z0) } 


exp{-j6t*r Jd8t  (4-2-3) 


The  first  integral  corresponds  to  the  Hertz  potential  of  a 
Hertzian  dipole  within  an  unbounded,  anisotropic  medium,  the 
second  integral  represents  the  effect  of  the  air-forest  'interface. 
These  integrals  were  evaluated  in  [42]:  the  first  integral 
describes  a  direct  wave  propagating  through  the  forest  to  the 
receiving  location;  the  second  integral  describes  a  wave  reflected 
from  the  air-forest  interface  and,  possibly,  a  lateral  wave 
propagating  along  the  air-forest  interface. 


‘valid  for  arbitrary  z  and  z  within  th«  slab. 
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Consider,  now,  the  Hertz  potential  within  the  anisotropic 
forest  slab  [Equation  (4-1-10)].  if  the  denominator  of  the 
integrand  is  expanded  in  a  power  series  using  the  relation 


(1  -  2)  -  1  +  z  +  z2  +  ... 


(  | z  J  <1 )  (4-2-4) 


Equation  (4-1-10)  may,  upon  interchanging  the  order  of  integra¬ 
tion  and  summation,  be  recast  into  the  form 


*11/ 

Az  (E'Z>  = 


-Id£ 


8TT2u)e0e 


where 


4 

£  s 

j=i  3 


(4-2-5) 


=  4  -  £  I« 


m=0 


and 


lm 


-/ 


“  ^  y  ^  x 


(4-2-6) 


*  exp{-j$t.r}d6t 


(4-2-7) 


*2m  =/C<St'»g<St) 


m, -  .  exP{-jT 2[2(m+l)H-(z+z.)]} 


•  exp{-jet*e)det 


(4-2-8) 


•  exp{-jBt*r)d8t 


(4-2-9) 


4m 


•-  Ar-MC'is.)  exp(~jT2i2(mti)H 


-|z-z.  I  ]  } 


•  exp(-j0t.r}d8t 


(4-2-10) 
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It  can  be  inferred  by  analogy  with  Equation  (4-2-3)  that  the 
integral 


‘10 


■/ 


exp{ - jx 0  | z-z9  |  } 


exp{-jBt*r}d3t 


(4-2-11) 


describes  a  direct  wave  propagating  through  the  forest  slab; 
the  integral 


20 


-fw 


exp{-jx. ( 2H-z-z  0) ) 


exp{-j£t *r}d^t  (4-2-12) 


describes  a  wave  reflected  from  the  air-forest  interface  and, 
possibly,  a  lateral  wave  propagating  above  the  air-forest  inter¬ 
face  as  well.  Further  analogy  with  Equation  (4-2-12)  suggests 
that  the  integral 

/exp{-jx, (z+z0) } 

Ra (-t  - — -  exp{-jBt.r}dBt  (4-2-13) 

describes  a  wave  reflected  from  the  forest-ground  interface; 
any  lateral  wave  that  might  be  propagated  below  this  interface, 
however,  would  be  highly  attenuated  by  ground  losses  and  so 
may  be  neglected. 


More  generally,  the  integrals  Ijm  can  be  evaluated  asymp¬ 
totically  using  the  method  of  steepest  descents,  just  as  Equation 
(4-2-12)  was  evaluated  in  [1,  Section  5.3}.  This  method  leads 
to  a  saddle-point  result  for  each  integral  as  well  as  a  branch- 
cut  contribution  if  a  branch-cut  is  crossed  by  the  steepest 
descent  path.  The  saddle-point  results  for  the  Ijm  integrals 
are 


I(r)  -  - 
lm 


a  lm  g  lm 


expf-jk./e^  Rlm} 
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(4-2-14) 
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2m  a 
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expf-jk./^  R2m> 
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(4-2-15) 


I<r)  -  - 
3m 
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exp{-jk.R/e^  R3mJ 
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(4-2-16) 


(4-2-17) 


Iir)  ~  “2ifRm+l(e/  )R*"+t  ( 0  ) 

4m  a  4m  g  4m 


exp{-jk.R/FT  R 


z  4  m 


aR 


4m 


where 


cos9  -  /FT ( 1  -  e  sin2©)*5 

R,(0)  =  - 5 - i - f 

cos9  +  /e“(l  -  e  sin20 ) * 


(4-2-18) 


R„(0)  = 


e  cos0  -  e  sin20) 

_  g  _ t  z 


egcos0  +  /F^(2f^ezsin20)  * 


(4-2-19) 


are  the  Fresnel  reflection  coefficients  at  the  air-forest  and 
ground-forest  interfaces,  respectively,  and 


0lm  =  Arctan{r/a [2mH+| z-z» | ] } 


(4-2-20) 


02m 

Arctantr/aUfm+DH-tz+z*)  ] } 

(4-2-21) 

03m  * 

Arct an  { r/a [ 2mH+ ( z+z  0 ) 1 } 

(4-2-22) 

9  4m 

Arctan{r/a [2 (m+1) H- | ?-z0 | ] } 

(4-2-23) 

Rlm  = 

(r2  +  a2  [2mH+|z-z0  | )  J*5 

(4-2-24) 

R2m  = 

{r2  +  a2  [2  (m+1)  H- (z+z.)  ]  }** 

(4-2-25) 

R3m  * 

{r2  +  a2  [2mH+(z+z„)  ] 

(4-2-26) 

R4m 

{r2  +  a2  (2  (m+1)  H- |  z-z„  |  ] 

(4-2-27) 

Branch-cut  contributions  to  the  Ijm  integrals  must  also 
be  included  whenever  the  0jm  exceed  the  so-called  critical  angle 

0  =  Arcsin(l//e~)  (4-2-28) 


The  branch-cut  contributions  to  the  I ^  integrals  are 


I(b)  ~  -4ire  mRm(0  ) 
lm  z  g  c 


exp{-jk0  [r+a(e  -l)ls[2mH+jz-z0  |]  ]  } 


k0(z-l)r: 


(4-2-29) 


exp{- jk„  [r+a (e  -1)  **  [2  (m+i)H-(z+z„)  ]  ]  } 


k0(ez-l)r2 


(4-2-30) 


~  -4Tie  mRm+1  (0  ) 
3m  z  g  c 


exp{-jk0 [r+a (ez~l) ^ [2mH+(z+z0) ] ] } 
k0 (e2-l)r 


(4-2-31) 


*lm  '  ->’S<m+1)Rg+,(0c) 


exp{- jk„  [r+a (e  -1)  **  [2  (m+1)  H-  |  z-z0  [  ]  ]  } 


k0  (e  -1) r 


(4-2-32) 


More  accurate  representations  for  the  branch 


■jm 


can  be  inferred  by  analogy  from  Equat 


-cut  contributions 
ions  (5-3-37)  and 


(5-3-40)  of  [42],  Constraints  bearing  upoif  the  validity  of 
these  results  can  also  be  inferred  from  [42], 

Because  branch-cut  Contributions  are  included  only  when 

9jm  >  9c  I  (4-2-33). 

the  number  of  branch-cut  contributions  is  finite.  As  a  conse¬ 
quence,  Equation  (4-2-6)  can  be  re-written  in  tne  form 


s.  =  T  i!r)  +  Y-  i! 

3  m=0  3in  m=0  3m 
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(b) 


The  electromagnetic  fields  (J2  and  II)  with! 
forest  slab  can  be  derived  from  the  Hertz  ve 


(4-2-34) 


in  the  anisotropic 
ztor  potential  A*1 
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[Equation  (4-2-5)]  by  using  Equations  (4-1-4)  and  (4-1-5).  For 
example,  the  vertically-polarized  electric  field  vector  E  is 


Ez  =  E^d)  +  E{zr)  +  E{zn 


(4-2-35) 


where 


^d)  (r,z)  =  sin20 


exp{-jk0/I^  R1Q} 
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(4-2-36) 
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(4-2-38) 


In  Equation  (4-2-35),  E^d)  represents  a  direct  wave  propagating 
through  the  anisotropic  forest  slab  from  the  transmitter  to 
the  receiver;  represents  multiply-reflected  waves  also 
propagating  through  the  anisotropic  forest  slab  from  the  trans¬ 
mitter  to  the  receiver  (refer  to  Figure  4-3a];  and  E^fc)  represents 
multiply-reflected  lateral  waves  associated  with  the  air-forest 
interface  (refer  to  Figure  4-3b] .  The  first  summation  in  Equation 
(4-2-37)  corresponds  to  waves  which  experience  an  even  number 
of  reflections  in  propagating  from  the  transmitter  to  the  receiver 
and  for  which  the  first  reflection  is  from  the  ground;  the  second 
summation  corresponds  to  an  odd  number  of  reflections ,  the  first 
of  which  is  from  the  air;  the  third  corresponds  to  an  odd  numjber 
of  reflections,  the  first  of  which  is  from  the  ground;  and  the 
fourth  summation  corresponds  to  an  even  number  of  reflections, 
the  first  of  which  is  from  the  air. 


4.3  Intra-Forest  Multipath 

The  anisotropic  forest  slab  model  differs  from  the 
anisotropic  forest  half-space  model  by  accounting  for  ground 
reflections.  The  principal  consequence  of  ground  reflections 
is  the  appearance  of  intra-forest  multipath.  Two  types  of 
intra-forest  multipath  can  be  identified:  multiply-reflected 
space  waves,  described  by  Equation  (4-2-37)  and  illustrated 
in  Figure  4-3 (a);  and  multiply-reflected  lateral  waves,  described 
by  Equation  (4-2-38 )and  illustrated  in  Figure  4— 3 ( b ) . 
Intra-forest  multipath  can  either  enhance  or  degrade  the  received 
signal  level  according  to  whether  the  concomitant  phasor 
interference  is  constructive  or  destructive;  in  addition, 
intra-forest  multipath  can  lead  to  pulse  distortion  and 
inter-symbol  interference. 


It  is  apparent  from  Equation  (4-2-37)  and  (4-2-38)  that 
the  Fresnel  reflection  coefficients  at  the  air-forest  and  forest- 
ground  interfaces  [denoted  by  Ra(0)  and  by  R  (e  )  ,  respectively] 
play  an  important  role  in  determining  the  relative  contributions 
of  the  contributing  multipath.  Thvi  behavior  of  these  reflection 
coefficients  is  examined  in  the  following  sub— section.  Successive 
sub-sections  address  the  intra-forest  multipath  associated  with 
multiply-reflected  space  waves  and  multiply-reflected  lateral 
waves . 


4*3.1  Fresnel  Reflection  Coefficients 


The  Fresnel  reflection  coefficients  associated  with  the 
air-forest  and  forest-ground  interfaces  are  given,  respectively, 
by  ' 

cosO  -  /F^(l  -  ezsin20)Jj 


R_  ( <) ) 
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cos0  +  /e  (1  -  e  sin20) 
t  z 


(4-3-1) 
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l<0>  ’  "a - =~F, - - - r 

E_cos0  +  /eJi1-  e  sin 2 0  ) ^ 


(4-3-2) 


g  “t  z 

where  efc  and  ez  are,  respectively,  the  transverse  and  longitudinal 
components  of  the  uniaxial  effective  dyadic  permittivity  of 


the  forest;,  eg  is  the  permittivity  of  the  ground;  and  0  is  the 
incidence  angle  measured  relative  to  the  normal. 


The 


magnitude  |  R  |  and  the  phase  |  Arg  |  of  the  Fresnel 


reflection  coefficient  for  the  fores -ground  interface  R  (0) 

g 


is  shown  in  Figure  4-4  (a)  as  a  function  of  the  glancing  angle 
[complement  of  the  incidence  angle]  for  the  case  of  moist  ground 


<eg  _30»  °g  =  °-05  s/m)  and  a  frequency  of  300  MHz.  At  the 
Brewster  angle,  |  R  |  exhibits  a  pronounced  minimum  and  the  phase 

i  _  i 


|Arg|  changes  by  nearly  180  degrees.  At  grazing  incidence 


(glancing  angle  equal  to  zero)  Rg(90°)  =  -1.  Analogous  results 


are  shown  in  Figure  4-4  (b)  for  the  case  of  dry  ground  (e1 
i5,  og  *  0.005  S/m).  Although  Equation  (4-3-2)  reveals  that 
that  Pg(0)  also  depends  upon  the  effective  dyadic  permittivity 


of  the  forest,  because  e  and  e  differ  from  unity  by  only  a 
few  tenths  of  one  percent,  Rg(0)  is  virtually  independent  of 
the  biophysical  forest  parameters. 

The  magnitude  |  R|  and  the  phase  | Argj  of  the  Fresnel  reflec¬ 
tion  coefficient  for  the  air-forest  interface  R  (0)  is  shown 
in  Figure  4-5(a)  for  the  case  of  lossless  leaves  (  =  40, 

e u"  =  0)  and  a  frequency  of  300  MHz.  Note  that  the  corresponding 

effective  dyadic  permittivity  is  purely  real  (  e."  =  e  "  =  0). 

z 

For  glancing  angles  less  than  the  critical  glancing  angle, 

0  =  Arcsin ( 1/ /e  )  (4-3-3) 

c  z 

| R |  is  unity  and  the  incident  radiowave  experiences,  total  internal 
reflection.  However,  whereas  at  grazing  incidence  (glancing 
angle  equal  to  zero)  R  (90°)  =  -1,  at  the  critical  glancing 

cl 

angle  R  (G  )  =  +1.  For  glancing  angles  greater  than  critical, 

cl  c 

( R |  is  nearly  zero.  At  the  Brewster  angle  |r|  exhibits  a  barely 
discernible  minimum  and  the  phase  |Arg|  changes  abruptly  by 
180  degrees.  Analogous  results  are  shown  in  Figure  4-5 (b)  for 
the  case  of  lossy  leaves  1  Note,  however,  that  as  a  consequence 
of  ohmic  losses  within  the  leaves,  the  incident  radiowave  no 
longer  experiences  total  internal  reflection  for  glancing  angles 
less  than  critical;  further,  the  phase  change  associated  with 
the  Brewster  angle  is  no  longer  abrupt. 


Multiply-Reflected  Space  Waves 


The  intra-forest  multipath  associated  with  multiply- 
reflected  space  waves  is  described  by  Equation  (4-2-37)  and 
illustrated  in  Figure  4-3(a).  It  is  apparent  from  Equation 
(4-2-37)  that  the  relative  strength  of  the  intra-forest  multipath 
can  be  estimated  by  considering  the  relative  magnitudes  of  the 
composite  forest  reflection  coefficients  [e.g.,  R  ,  R  ,  R  R  , 

cl  3  ^ 

R  R.  R  R  R  ,  RPR  (refer  to  Figure  4— 3(a))]  "  gnitudes 

gaaga  gag 

of  the  composite  Fresnel  reflection  cot.  -vines  at  300  MHz 

for  a  forest  of  leaves  (FV  =  0.1%)  above  a  moist  earth  (e  '  * 

9 

30,  0^  =*  0.05  S/m)  have  been  plotted  in  Figure  4-6  as  a  function 


20  30  40  50  60  70  80  90 

Glancing  Angle  (degrees) 

-5(b):  Reflection  Coefficients  (Lossy  Leaves) 
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Figure  4-6 (b):  Reflection  Coefficients  (Composite)  dB 


of  path  length  using  the  expression 


0  = 

where,  H  = 

Zo  = 

z  = 

r  = 

e  = 

n  = 

to  relate  the  path  length  to  the  incidence  angle.  It  is  apparent 

from  this  figure  that  for  these  biophysical  forest  parameters 

and  a  frequency  of  300  MHz  the  ray  reflected  only  once  from 

the  air-forest  interface  Ithe  only  ray  accounted  for  in  the 

anisotropic  half-space  model)  is  important  at  all  distanqes 

in  excess  of  about  300  itieters.  The  ray  reflected  once  from 

the  ground  is  important  not  only  at  ,  distances  in  excess  of  300 

meters,  bur  ;lso  at  very  short  distances  as  well.  The  R  R 

a  g 

and  RgRa  contributions  become  important  only  at  ranges  exceeding 

about  300  meters;  and  the  R  R  R  and  R  R  R  contributions  can 

a  g  a  gag 

probably  be  neglected  for  ranges  less  than  1000  meters.  It  should 
be  noted  that  Figure  4-6  does  not  account  for  antenna  pattern 
discrimination,  differential  geometric  spreading  losses,  or 
differential  absorptive  losses. 

4.3.3  Multiply-Reflected  Lateral  Waves 

The  intra-forest  multipath  associated  with  multiply- 
reflected  lateral  waves  is  described  by  Equation  (4-2-38)  and 
illustrated  in  Figure  4-3<b).  The  correspondence  between  the 
multiply-reflected  lateral  waves  of  Figure  1-3(b)  and  the  summa¬ 
tions  of  Equations  (4-2-38)  is  summarized  in  Table  4-1.  The 

•The  upper  (lower)  quantity  within  j  J  is  to  be  used  If  tne  radiowave 
reflection  nearest  the  transmitter  (receiver]  is  from  the  air  (ground). 


Arctan 


(n-l)H  +  1h-z o I  +  |H~  z  |  j 
I  Zo  I T  j  z  !  R 


(4-3-4) 


height  of  forest 

height  of  transmitter  above  forest  floor 
height  of  receiver  above  forest  floor 
distance  between  transmitter  and  receiver 
incidence  angle  of  reflected  radiowave 
number  of  reflections  experienced  by  radiowave 


Table  4-1:  Multiply- Ref lected  Lateral  Wave  Ray  Paths 
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multipath  geometry  illustrated  in  Figure  4-3 (b)  can  be  easily 
constructed  by  locating  the  appropriate  transmitter/receiver 
images  arising  witn  successive  radiowave  reflections  from  the 
air  and  ground  interfaces  at  the  critical  angle  0c-  The  location 
of  these  images  relative  to  the  air-forest  interface  can  be 
inferred  from  the  exponents  appearing  in  the  summations  of  Equa¬ 
tion  (4-2-38);  those  corresponding  to  Table  4-1  are  shown  in 
Table  4-2. 


Equation  (4-2-33)  can  be  used  in  consonance  with  Equations 
(4-2-20)  through  (4-2-27)  to  show  that  associated  with  each 
(j,m)  component  lateral  wave  is  a  minimum  excitation  distance 


r  . 
min 


(4-3-5) 


where  2jm'  the  so-called  [65]  separation  distance  is  given  by 

zlm  =  2mH  +  lz-'“l  (4-3-6) 
z2m  =  “  (z+2<>)  (4-3-7) 
z3in  =  2mH  +  (7!-Zo)  (4-3-8) 
z4m  =  2(m+l)H  -  |z-z0|  (4-3-9) 


The  separation  distance  [z.  ],  the  minimum  excitation,  distance 

3m 

[rminl»  and  the  relative  strength  of  the  multiply- ref lected 
lateral  waves  have  been  calculated  for  several  1000-metex ,  300-MHz 
radio  links  situated  within  a  forest  of  leaves  (FV  =  0.01%) 
and  above  a  moist  ground  (e^  =  30,  a  =  0.05  S/m)  and  presented 
in  Table  4-3.  Note  that  the  non-ground-reflected  lateral  wave 
[the  only  lateral  wave  accounted  for  in  the  anisotropic  half-space 
model)  is  dominant.  The  phasor  sum  of  all  contributing  lateral 
wave  components  has  been  calculated  and  its  relative  strength 
denoted  in  the  Table  by  LDBT.  As  expected,  it  differs  only 
slightly  from  that  of  the  non-ground-reflected  lateral  wave. 
When  the  number  of  contributing  lateral  waves  is  large,  the 
phasor  sum  is  well  approximated  by 


Table  4-3 


LATERAL  WAVE  MULTIPATH 

PATH:  Forest  Ht.  =  25. m  Ht-xmtr  =  22.5m  Ht-rcvr  =  22.5m 

Critical  Angle  =  3.83  Deg.  Reflection  Coef  =  .4644 
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Excitation 

Distances 

(m) 
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5.00 
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Separation  Distances  (m) 
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Relative  Field  Strength 
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-16.58 
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.00 

-11.66 
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-16.58 
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LDBT  = 

-.16338  dB  LDBTT  *  -4 

.32773  dB 

K  (  z  ,  z  o  ,  H  ) 


(4-3-10) 


[1  +  R  (9  )  exp{  jicz  o  }  ]  [  1  +  R  (0  )exp{j«rz}] 

_ q  c  _  q  c 

[1  +  Rg  (  ©c  )exp{  j k II >  ]  2 

where 

k  =  2k0/a/fJ  “IT  (4-3-11) 

z 

This  too  has  been  calculated  and  is  denoted  by  LDBTT  in  the 
Table.  It  is  apparent  from  this  calculation  that  when  only 
a  few  lateral  wave  components  contribute  significantly,  Equation 
(4-3-10)  is  inappropriate. 

4 . 4  Transmission  Loss 

The  most  important  parameter  characterizing  radiowave 
propagation  channels  is  transmission  loss.  Loosely  defined, 
transmission  loss  is  the  ratio  of  received  power  to  transmitted 
power  (or  its  inverse).  Basic  transmission  loss,  however,  is 
specifically  defined  (for  unmodulated  r-f  carriers]  as  the  ratio 
of  the  power  transmitted  to  the  power  received  as  measured  at 
the  terminals  of  isotropic,  lossless,  co-polarized  antennas. 
The  concept  of  basic  transmission  loss  may  be  extended  to 
wideband,  modulated  r-f  carriers  by  appropriately  weighting 
the  basic  transmission  loss  by  the  transmitted  power  spectral 
density. 

4.4.1  Basic  Transmission  Loss 

For  radiowaves  propagating  through  freespace,  the  basic 
transmission  loss  (in  decibels)  can  be  calculated  from  the  well- 
known  expression 

Lfg  =  32.4  +  2  0  £og  i  o  dj^jYi  +  20£ogjofMHz  (4  —  4  —  1) 

where  dkm  is  the  distance  between  the  antennas  in  kilometers 
and  f mhz  the  radiowave  frequency  in  Megahertz.  This  equation 
has  been  used  to  plot  in  Figure  4-7  the  frequency  dependence 
of  basic  transmission  loss  in  free  space  over  a  1000-meter  path. 
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In  calculating  the  basic  transmission  loss  associated  with 
the  stratified  forest  models  described  previously,  it  is  important 
to  recognize  that  postulated  for  those  models  was  a 
vertically-polarized  electric  dipole  antenna.  Such  an  antenna 
is  not  isotropic.  Nevertheless,  the  basic  transmission  loss 
can  still  be  determined  by  noting  that:  (1)  in  free-space, 

the  vertically-polarized  electric  field  can  be  determined  from 
Equation  (4-2-36)  by  setting  =  ez  =  1  so  that 

Efs  =  (u)Uo/4tt  )Id£  sin29!  0---P^k°R]'a-  (4-4-2) 

*1  0 

and  (2)  the  power  received  is  directly  proportional  to  the  field 
intensity 

I(x)  =  E(x) -E*(x)  (4-4-3) 

Thus,  for  the  stratified  forest  models,  the  basic  transmission 
loss  can  be  expressed  by 

Lb ( d , r , £ )  —  Lfg  +  20£og ; o | Ef g/Ed r r, Z I  (4—4—4) 

where  E d,r,l  is  defined  in  Equations  (4-2-36),  (4-2-37),  and 

(4-2-38).  Equation  (4-4-4)  has  been  used  in  Section  5.2  to 
assess  the  sensitivity  of  the  basic  transmission  loss  to  the 
biophysical  forest  parameters. 

4.4.2  Wideband  Basic  Transmission  Loss 

Consider  a  modulated  radiowave  having  a  transmitted  power 
spectral  density  S(w).  The  total  power  transmitted  is  then 

Prp  =  J's  (u>)du>  (4-4-5) 

If  the  basic  transmission  loss  associated  with  an  unmodulated 
r-f  carrier  of  angular  frequency  u  is  denoted  by  L(u>)  then  the 
total  power  received  is 

Pr  —  /s  ( w )  L  ( cu  )dw 


(4-4-6) 


* 


The  wideband  basic  transmission  loss  can  then  be  defined  as 


^wb  “ 


Pr  y*S(u)  )L(w  )do 

Pm  ~T  S  ( a)  )da) 


(4-4-7) 


or,  expressed*  in  decibels. 


Ljg  —  -lOfcog i  olvb 


(4-4-8 ) 


4  -  5  Doppler  Frequency  Shift 

The  stratified  forest  model  for  the  coherent  field  components 
of  the  propagating  radiowave  has  been  extended  to  account  for 
Doppler  shift  induced  by  terminal  motion.  This  development 
is  described  below.  Doppler  shift  of  the  coherent  component 
induced  by  motion  of  the  forest  biomass  (primarily  leaves  and 
branches)  is  anticipated  to  be  relatively  small  by  comparison. 

The  effects  of  terminal  motion  and  the  Doppler  shift  induced 
thereby  can  be  addressed  most  simply  by  first  considering  only 
the  direct  field  component 

,  exp{-jk0  /ez  Rd} 

Ez(d)  =  /a  — 2)  IdJt*sin20d  — - — — - - — =—  (4-5-1) 

\  4"/  a  Rd 


te) 


ld£.*sin20d 


(4-5-1) 


As  a  consequence  of  terminal  motion,  Rd  (the  distance  between 
the  transmitter  and  receiver)'  becomes  time-variant.  If  the 
terminal  motion  is  uniform  (no  acceleration),  Rd  may  be 

represented  by  the  first  two  terms  of  its  Taylor  series  expansion 

. 

so  that 


Ra  -  ?d+  ^ 


(4-5-2) 


Over  relatively  short  time  intervals  where 


6t  << 


[dRd/dt ] 


=  Atr 


(4-5-3) 


*  The  negative  sign  has  been  inserted  for  consistency  with  Equation  (4-4-4), 
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Equation  (4-5-1)  can  be  well  approximated  by 

/ coyri\  exp(  j<  fz~  (R°  +  R't)} 

Ez(d)  =  /a  1— o  z  d  d 


(£) 


Id  £  •  sin2  Sjj  • 


(4-5-4) 


Remembering  that  the  instantaneous  frequency  of  a  waveform 
proportional  to  exp{$(t)}  is  defined  as 


f  _  1  d<T 
2?  dt 


(4-5-6) 


where  <p"  represents  the  imaginary  part  of  <p,  and  remembering, 
too,  that  the  factor  exp{j2irf0t}  is  assumed  but  suppressed  in 
Equations  (4-5-1)  and  (4-5-4),  the  instantaneous  Doppler  shift 
in  frequency  of  the  direct  field  component  Ez(d)  relative  to 
the  carrier  frequency  fQ  is  given  by 


Since 


Afd  =  - T0Re{/ez}R^/c 


Rd  =  [r2  +  a2 (z  -  z0) 2 ] 


(4-5-6) 


(4-5-7) 


then 


9  (z  -  zQ) 


3(z  -  zD) 
3t 


(4-5-8) 


Rd  =  vrpined  +  Vzcos0d  (4r-5-9) 

where, 

vr  ~  radial  component  of  receiving  terminal  velocity 
relative  to  transmitting  terminal  velocity 

[=  9r/9t ] . 

V z  =  vertical  component  of  receiving  terminal  velocity 
relative  to  transmitting  terminal  velocity 

[=  9 ( z  -  zQ )/9t ] .  y 

Note  from  Equation  (4-5-6)  and  (4-5-9)  that  the  Doppler  shift 
is  directly  proportional  to  the  carrier  frequency  fQ  and  to 
that  component  of  the  differential  terminal  velocity  directed 
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between  the  terminals.  The  negative  sign  in  Equation  (4-5-6) 
signifies  that  relative  terminal  motion  tending  to  increase 
Rd  will  shift  the  received  frequency  downward  relative  to  the 
carrier. 

In  a  similar  manner,  the  instantaneous  Doppler  shift  of 
any  reflected  field  component  Ez(r)  relative  to  the  carrier 
frequency  fQ  can  be  found  to  be 


.Af  =  -f  Re{/e  }R'/c 
r  o  z  r 


(4-5-10) 


where 

R'r  =  Vrsin0r  -  Vzcos9r  (4-5-11) 


and,  in  addition  to  the  previously  defined  variables, 

Vz  ='  vertical  component  of  receiving  terminal  velocity 
relative  to  the  transmitting  terminal  image 
velocity  [=  3  ( z-rz0  )/3t  ] . 


The  instantaneous  Doppler  shift  of  any  lateral  wave  field 
component  E z(£)  relative  to  the  carrier  frequency  fQ  is  given 
by 


Af  £  »  -(f0/c)’{Vr  -  Re{a(ez  -  1 )  *5 }  -  V2  ]  (4-5-12) 


The  order  of  magnitude  of  the  Doppler  shift  may  be  estimated 
by  considering  two  vehicles  moving  apart  along  the  forest  floor 
at  speeds  of  100  kilometers  per  hour  and  utilizing  a  carrier 
frequency  of  600  megahertz.  The  Doppler  shift  of  the  direct 
field  component  will  dominate  and  will  be  approximately  equal 
to 

fd  =  -(600x10s) 

3600 

=  -100  Hertz. 

For  a  broadband  communication  system  utilizing  bandwidths  in 
excess  of  100MH  ,  Doppler  shifts  of  this  magnitude  can  be  expected 
to  have  negligible  effect  on  performance. 


^/(3xl03 


(4-5-13) 
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5 . 0  Sensitivity  of  the  Stratified  Forest  Model 


Under  this  task,  the  stratified  forest  model  has,  been  used 
to  investigate  the  sensitivity  of  the  basic  transmission  loss, 
the  (mean)  forest  pulse  response,  and  the  differential  pulse 
delay  to  variations  in  the  salient  biophysical  parameters  of 
the  forest  [refer  to  Section  2]  and  to  the  frequency, 
polarization,  and  path  length  of  the  radio  link.  Three  forest 
types  have  been  considered  to  illustrate  the  sensitivity  of 
the  stratified  forest  model:  an  uneven-aged,  trunk-dominated 
forest  characterized  by  a  trunk  number  density  (pt)  of  1000 
trunks  per  hectare  and  by  an  exponential  trunk-diameter 
probability  density  function*  having  a  mean  trunk  diameter  of 
6.35  centimeters  (2.5  inches);  an  all-leaf  forest  characterized 
by  a  fractional  volume  (FV)  of  0.1  percent  [5  centimeter  leaf 
radius,  1  millimeter  leaf  thickness,  and  133  leaves  per  cubic 
meter];  and  a  mixed-forest  of  trunks  and  leaves  derived  from 
the  other  two.  Branches  have  not  yet  been  incorporated  into 
the  stratified  forest  model. 

Although  all-leaf  forests  will  be  studied  systematically 
for  later  comparison  with  uneven-aged,  trunk-dominated  forests, 
it  seems  advisable  to  start  this  sensitivity  analysis  of  the 
stratified  forest  model  with  the  one  all-leaf  forest  considered 
previously  in  Reference  [42  ] .  This  will  permit  those  previous 
results  to  be  discussed  in  a  manner  not  possible  at  that  time 
and  to  support  an  evolutionary  development  of  the  formats  to 
be  used  in  the  systematic  sequel.  Figures  6-2  and  6-3  of  Re¬ 
ference  [42]  have  been  reproduced  here  in  a  newer  and  more 
complete  form  as  Figures  5-1  and  5-2  As  before,  the  number 


‘Although,  as  noted  in  Section  2.2.1,  the  upper  and  lower  limits  on  the  trunk 
diameters  .  of  real  forests  must,  in  practice,  be  reconciled  with  the  seni- 
infinite  limits  available  to  the  exponential  probability  density  function, 
for  present  purposes  the  exponential  probability  density  function  was  defined 
only  over  a  1  tol2  inch  diameter  range  at  discrete  1-inch  intervals.  The 
area  under  the  density  function  was,  of  course,  normalized  unity. 
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density  of  the  leaves  has  been  taken  as  133*  leaves  per  cubic 
.ueter  (corresponding  to  a  fractional  volume  of  0.1  percent]. 
The  transmitter  and  receiver  antennas  have  been  situated  3  5+ 
meters  below  the  forest  top  and  the  range  has  been  taken  as 
10C0  meters.  The  forest  pulse  response  shown  in  Figure  5-2 
has  been  determined  for  a  5-nanosecond  rectangular  pulse  having 
an  r-f  carrier  frequency  of  600  MHz. 

The  forest  transfer  functions  shown  in  Figure  5-1  have 

been  normalized  to  0  dB  for  the  direct  wave  at  600  MHz.  The 
direct  wave  is  reduced  by  interference,  with  the  reflected  wave 
at  all  frequencies,  with  nearly  complete  cancellation  occurring 
near  400  MHz;  their  resultant  is  less  than  the  lateral  wave 
at  all  frequencies.  Thus,  the  forest  transfer  function  of  the 
total  wave  follows  roughly  that  of  the  lateral  wave  with 
interference  lobes  arising  as  a  consequence  of  further 

interference  between  the  lateral  wave  and  the  direct-plus- 
reflected  wave.  The  forest  pulse  response  of  Figure  5-2  shows 

the  same  partial  cancellation.  The  lateral  wave  arrives  earlier 
than  the  direct  wave  because  it  travels  faster  in  free  space 

above  the  forest  than  the  other  waves  travel  through  the  forest. 
The  corresponding  pulse  envelopes  are  shown  in  Figure  5-3. 

The  subsequent  presentation  of  the  sensitivity  analysis 

involves  two  significant  changes  in  format  as  well  as  variations 
in  forest  type  and  path  geometry.  The  first  is  the  inclusion 
of  basic  transmission  loss.  This  is  discussed  in  Section  4-4 

and  in  Section  5-1  below.  The  second  is  the  insertion  of  an 
ideal  (rectangular)  bandpass  filter  centered  on  the  carrier 

frequency  of  600  MHz  to  improve  the  simulation  of  the  radio 

equipment.  For  an  ideal  rectangular  pulse  with  its  sine-shaped 


*Al 1  curves  presented  in  Reference  (421  citing  a  leaf  number  density  of  200 
leaves  per  cubic  meter  were,  in  fact,  based  upon  a  leaf  number  density  of 
(2/3) (200)  =  133  leaves  per  cubic  meter. 

+The  radio  path  geometry  is  slightly  different  frem  that  actually  used  in 
Reference  (42).  There,  the  transmitter  and  receiver  antennas  were  situated 
3  and  4  meters,  respectively,  below  the  forest  cop  (not  4  and  6  meters, 
respectively,  as  stated);  here,  both  transmitter  and  receiver  antennas  are 
situated  3.5  meters  below  the  forest  top. 
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Figure  5-1 s  Transfer  Function  for  Leaf  Forest  (R  =  1000m,  S/2 
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Figure  5-3:  Pulse  Envelope  for  Leaf  Forest 
(R  -  1000m,  S/2  -  3.5m) 
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amplitude  spectrum,  the  high-  and  low-frequency  spectral  skirts 
led  to  intolerable  pulse  distortion.  Therefore,  in  the  sequel 
only  the  central  lobe  of  the  pulse  spectrum  is  passed  by/ideal 
filter.  The  corresponding  transmitted  pulse  shape,  no  longer 
rectangular,  is  shown  in  Figure  5-4. 


5.1  Basic  Transmission  Loss 

The  basic  transmission  loss  experienced  within  an  all-leaf 
forest  characterised  by  a  fractional  volume  (FV)  of  0.1  percent 
is  presented  in  Figures  5-5  and  5A-1  to  5A-8.  The  vertically-polarized 
transmitter  and  receiver  antennas  are  situated  2.5  meters  below 
the  forest  top  and  separated,  successively,  by  200,  500,  and 

1000  meters.*  The  basic  transmission  losses  associated  with  the 
direct  and  reflected  waves  are  virtually  indistinguishable  and 
for  frequencies  less  than  about  500  MHz  increase  about  6  dB/octave 
of  frequency.  The  direct  and  reflected  waves  destructively 
interfere  so  that  their  resultant  is  nearly  20  decibels  below 
either  of  them  at  R  ■  1000  m.  Particularly  strong  destructive 
interference  is  apparent  near  500  MHz.  The  basic  transmission 
loss  of  the  total  at  1000  m  is  clearly  dominated  by  that  of 
the  lateral  wave  which  decreases  about  12  dB/octave  of  frequency. 
Note  that,  in  general,  a  leaf-dominated  forest  behaves  as  a 
low-pass  filter  for  vertically-polarized  radiowaves. 

The  losses  for  the  same  ranges  with  S/2  ,  *  5  m  are  shown 
in  Fiqures  5A-4  to  5A-6,  and,  for  S/2  ?  10  n  in  Figures  5A-7  to 
5A-9.  The  chief  points  to  notice  are  the  absence  of  the  lateral 
waves  as  the  angle  of  incidence  exceeds  the  critical  angle  (see 
Section  4.2.6). 

The  basic  transmission  loss.,  for  the  trunk-dominated  forest 
is  shown  in  Figures  5A-10  to  5A-18.  As  for  the  leaf-dominated  forest, 
the  basic  transmission  losses  associated  with  the  direct  and  re¬ 
flected  waves  are  virtually  indistinguishable  and  they  destruc¬ 
tively  interfere  so  that  their  resultant  lies  10-15  d9  below 
them.  The  basic  transmission  loss  of  the  total  is  clearly 
dominated  by  that  of  the  lateral  wave  (they  are  indistinguishable) 

which  increases  approximately  6  dB/octave  of  frequency. 

\  .  ■ 

•except  for  Figure  5-5,  for  which  S/2  ■  3.5ra  and  R  ■  1000m  as  in  Figure  5-1. 
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5—4:  Filtered  Pulse  Shape 
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5 . 2  Forest  Pulse  Response 


The  mean  forest  pulse  response  has  been  determined  fob 
a  wideband  digital  radio  system  having  a  600  MHz  carrier  frequency 
and  a  300  MHz  r-f  bandwidth  as  above.  The  pseudo-noise  chips 
are  derived  by  passing  5-nanosecond  rectangular  pulses  through 
an  ideal  bandpass  filter.  The  foresc  pulse  response  for  the 
leaf-dominated  forest  is  shown  in  Figures  5B-1  to  5B-3  for 
S/2  =  2.5  m,  5B-4  to  5B-6  for  S/2  =  5  m  and  5B-7  to  5B-9  for 
S/2  =  10  m.  The  forest  pulse  responses  associated  with  radiowave 
propagation  along  the  direct  path  are  shown  uppermost  in  these 
figures.  The  propagation  delay  along  the  direct  path  x^,  observed 
to  be  about  3340  nanoseconds  for  R  =  1000  m,  can  also  be  estimated 
from  the  approximation 


Td 


c 


( 5-2-1. ) 


where  t'z  is  the  real  part  of  the  longitudinal  component  of  the 
effective  dyadic  susceptibility  [refer  to  Figure  3-1 (b) ] .  For 
this  example,  =  1000  m,  z'z  =  1.004265,  and  c  =  3x10®  m/sec. 
The  companion  pulse  response  for  the  reflected  pulse  also  shown 
in  Figure  5B-3  exhibits  essentially  the  same  delay.  The  rather 
shallow  glancing  angle  at  the  air-forest  interface  associated 
with  this  particular  geometric  configuration  (about  0.14  degrees) 
ensures  a  relatively  strong  reflected  pulse  (at  600  MHz  the  reflec¬ 
tion  coefficient  is  0.9942),  although  essentially  in  antiphase 
with  the  direct  pulse.  As  shown  in  the  figure,  the  direct  and 
reflected  pulses  virtually  cancel  each  other  and,  as  a 
consequence,  the  contribution  of  the  lateral  wave  dominates 
the  resultant.  The  propagation  delay  associated  with  the  lateral 
wave  is  essentially  that  of  free  space. 

The  forest  pulse  responses  for  the  trunk-dominated  forest  are 
shown  in  Figures  5B-10  to  5B-18.  In  Figure  5B-12,  the  pulse  response 
associated  with  the  direct  wave  arrives  prior  to  the  pulse 
response  associated  with  the  lateral  wave.  This  is  a  consequence 
of  the  fact  that  the  longitudinal  component  of  the  effective 
dyadic  susceptibility  is  negative.  Equation  (5-2-1)  can  still 


be  used  to  estimate  the  arrival  time  of  the  direct  pulse,  although 
for  this  example,  because  z'z  =  0.99903  the  velocity  of  pulse 

propagation  appears  to  exceed  that  of  light  in  vacuum.  A  more 
exact  expression  for  pulse  delay  time  and  the  Kramers-Konig 
relation  can  be  employed  to  refute  this  contention  in  general; 
the  apparent  paradox  that  Figure  5B-12  supports  this  contention 
can  be  resolved  by  noting  that  the  strictly  band-limited  trans¬ 
mitted  pulse  cannot  be  localized  in  time.  The  small,  rapid 
oscillations  may  be  the  so-called  Sommerfeld  precursor-  [36]. 
in  any  case,  the  high  specific  attenuation  associated  with  the 
propagation  of  a  vertically-polarized  wave  through  a  trunk- 
dominated  forest  [refer  to  Figure  3-2 (a)]  severely  attenuates 
all  but- the  lateral  wave  which,  as  is  apparent  from  Figure  5B-12 
dominates  the  resultant  pulse  response. 
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Figure  5 A- 4;  Basic  Transmission  Loss  (Leaf  Forest) 
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Figure  5A-7:  Basic  Transmission  Loss  (Leaf  Forest) 
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Figure  5A-10:  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest) 


Figure  5A-11:  _  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest) 


Figure  5A-12:  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest) 


Figure  5A-13:  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest) 


Figure  5A-14:  Basic  Transmission  Lo’ss  (Inhomogeneous  Trunk  Forest) 


Figure  5A-15:  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest) 
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Figure  5A-16:  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest.) 
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Figure  5A-17:  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest) 


Figure  5A-18:  Basic  Transmission  Loss  (Inhomogeneous  Trunk  Forest) 
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Figure  5B-2:  Pulse  Response  (Leaf  Forest) 
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Figure  5B-7:  Pulse  Response  (Leaf  Forest) 
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Figure  5B-17:  Pulse  Response  {Inhomogeneous  Trunk  Forest) 
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6 . 0  Incoherent  Forest  Scatterin 

The  electromagnetic  fields  of  radiowaves  propagating  through 
a  forest  can  be  decomposed  int-o  mean  (coherent)  and  diffuse 
(incoherent)  components.  At  low  frequencies  (HF  and  below) 
the  diffuse  component  is  relatively  small  and  only  mean  fields 
need  be  considered.  As  the  frequency  becomes  higher,  spatial 
fluctuations  in  the  field  become  more  important,  and  the  diffuse 
component  of  the  field  must  be  taken  into  account.  In  this 
section,  the  behavior  of  the  diffuse  component  at  VHP  and  OHF 
frequencies  is  examined  for  a  trunk-dominated  forest. 

The  trunk-dominated  forest  consists  solely  of  tree  trunks 
which  are  parallel  to  each  other  but  randomly  placed  perpendicular 
to  the  forest  floor.  It  is  a  useful  model  since  trunks  are 
the  largest  forest  component  and  the  first  to  give  rise  to  sub¬ 
stantial  random  or  diffuse  component  of  the  fields  as  the 
frequency  is  increased.  In  addition,  if  the  trunks  are  assumed 
to  be  circular,  and  propagation  is  parallel  to  the  forest  floor. 
Maxwell's  equations  can  be  replaced  by  a  scalar  wave  equation; 
this  represents  a  substantial  simplification  in  the  mathematical 
analysis  required. 

To  characterize  propagation  in  the  forest  when  the  diffuse 
component  of  the  field  is  significant,  the  space-frequency  corre¬ 
lation  functicln  of  the  field  is  required.  This  is  the  correlation 
function  between  the  field  component  at  space  point  x  and  fre¬ 
quency  <*•  with  the  conjugate  of  the  field  component  at  space  point 
x  and  frequency  u.  An  approximate  equation  for  this  correlation 
function  can  be  obtained  for  the  two-dirtiensional  trunk  forest. 
This  equation  is  the  starting  point  for  most  of  the  calculations 
of  this  section. 

Once  the  correlation  equation  has  been  obtained  and  the 
scattering  properties  of  the  trunks  discussed  two  problems  are 
addressed.  The  first  is  the  solution  of  the  correlation  equation 
for  trunks  thin  compared  to  a  wavelength  with  x  ■  x  but  «  f  u. 
With  this  frequency  correlation  function  the  coherence  bandwidth 
of  the  channel  can  be  determined.  Following  this  calculation. 


correlation  function 


the  intensity  and  space  correlation  function  are  evaluated 
for  the  general  case  of  a  forest  of  resonant  trunks,  i.e.,  trunks 
that  are  not  necessarily  small  compared  to  a  wavelength.  The 
intensity  calculation  represents  a  generalization  of  that  for 
thin  trunks  that  appears  in  [  4 2 1  -  The  results  of  both 
calculations  yield  essential  physical  information  on  channel 
behavior  which  will  be  discussed  subsequently. 


6 . 1  Mean  and  Correlation  Equations 

In  this  section,  the  equation  for  the  mean  field  and  the 
equation  for  the  space-  frequency  correlation  function  are 

developed.  Before  this  can  be  done,  the  two-dimensional  forest 
problem  must  be  formulated.  Although  some  of  the  material  in 
this  section  has  already  appeared  in  [42],  it  is  repeated  here 
for  completeness  and  convenience. 

6.1.1  Formulation  of  Two-Dimensional  Trunk  Problem 

Consider  a  trunk-dominated  forest  represented  by  a  col¬ 
lection  of  parallel  circular  dielectric  cylinders  having  radius 
a  and  complex  relative  dielectric  constant  .  The  collection 
of  trunks  is  confined  to  stand  within  an  area  S  as  is  shown 

in  Figure  6.1.  If  only  an  electric  current  source  having  no 
variation  in  the.  direction  parallel  to  the  cylinders  is 

considered,  the  electric  field  obeys  tht*  following  scalar  wave 
equation 

[V£  +  k|c  (xt,w)  ]Ez(xt,w)  =  ju)U0Jz  (xt,io)  (6-1-1) 

Here  Ez(xt,u>)  is  the  only  component  of  the  electric  field  excited; 
it  depends  only  on  the  transverse  position  vector  x.t  35  xx’  +  yy°. 
In  the  above  formulation,  the  z  coordinate  is  parallel  to  the 

cylinders. 

The  current  source  can  within  the  forest  region 

S  or  outside  of  it.  In  some  cases  the  source  is.  allowed  to 
recede  to  infinity  in  such  a  way  that  it  generates  plane  waves 
normally  incident  on  the  forest.  The  two-dimensional  source 
distribution  is  not  as  realistic  as  a  threfe-dimensional  point 


source,  but  the  scalar  formulation  of  Equation  (6-1-1;  lead  * 
to  a  much  more  tractable  problem.  Since  a  point  source  can 
be  represented  as  a  superposition  of  phased  line  sources,  the 
results  can  be  generalized  to  the  three-dimensional  source  case 
by  employing  superposition.  This,  however,  will  not  be  done 
at  this  time. 

The  dielectric  characteristics  of  the  forest  are  specified 
by  the  relative  dielectric  permittivity,  e(xt,io).  This 
permittivity  can  be  expressed  as  a  sum  of  the  susceptibilities 
of  the  individual  scatterers  as  follows: 


IT 

[x  ,u>)  =  1  +  ]jT  x .  (x  ,< 
c  j=l  J  c 


(6-1-2) 


where  Xj  is  the  susceptibility  of.  the  jth  scatterer.  It  is 
convenient  to  express  the  Xj  as  a  translation  of  the  suscepti¬ 
bility,  X  ,  of  an  identical  scatterer  located  at  the  origin. 
If  Xtj  is  the  vector  from  the  origin  to  the  center  of  the 
scatterer  then 


where 


Xj(xt,tu)  =  X(xt  -  Xt_.,w) 


X(xt,w)  = 


x4(«) 


(6-1-3) 


|xfc|  -a 


I  x  t  i  >  a 


Here  X^'io)  is  the  susceptibility  of  the  bulk  tree-  trunk  material. 
It  is  assumed  to  be  isotropic. 

Following  the  methodology  of  Section  2.1.1,  142],  Equation 
(6-1-1)  is  expressed  in  operator  notation.  If  the  following 
quantities  are  defined:' 


L  =  -(V*  +  k|) 


vj  *  k5xj 


(6-1-4 ) 


y  =  e_ 


g  =  -ju>u„j_ 


(6-1-5) 


then  Equation  (6-1-1)  can  be  written  as 


(L  '  £  V:)T  =  9  ,6-1-6' 

This  notation  has  been  used  in  Section  2.2.3,  (42]  to  derive 

the  mean  equation.  It  will  be  used  in  this  report  to  derive 
the  space-frequency  correlation  equation  (Appendix  A).  It 
provides  a  concise  mechanism  for  writing  defining  equations 
for  the  mean  and  correlation  in  the  main  body  of  the  text. 

The  formulation  of  the  scattering  problem  thus  far  is  in 
terms  of  the  susceptibility  of  individual  scatterers  translated 
1  from  the  origin.  Past  experience  has  shown  that  it  is  more 
useful  to  express  the  characteristics  of  a  scatterer  by  its 
response  to  an  incident  field.  For  this  purpose,  the  transition 

operator  T  is  introduced.  If  is  the  field  incident  on  a 

scatterer  having  susceptibility  X(xt»w)  which  is  located  at 
the  origin,  then  the  induced  sources  geq  generated  within  the 
scatterer  are 

geq  =  T^i  =  -fti*,*'  ,«*>)^i(x’  ,oj)dx'  (6-1-7) 

Here  the  transition  operator  has  been  represented  in  terms  of 
an  integral  operator  having  kernel  t(x,x',u>).  From  a  knowledge 
of  t(3c,x.' ,“)  all  properties  of  the  scatterer  can  be  determined. 

It  has  been  shown  in  Section  2*2.2,  [42]  that  the  transition 

operator,  T(Xtj),  for  scatterers  located  at  Xtj  can  simply 

expressed  in  terms  of  the  transition  operator  of  the  scatter 
located  at  the  origin.  This  relationship  is  given  by 

T(Xtj)  ^  =  y*t  (x  -  Xtj  ,x'  -  Xtj  ,o»)  ,u)dx'  (6-1-8) 

Using  Equation  (6-1-8)  the  properties  of  all  particles  can  be 
expressed  solely  in  terms  of  t(x,x* ,  w) . 


The  scatterer  can  be  characterized  in  an  alternative  manner 
by  specifying  it's  far-field  response  to  a  unit  amplitude  incident 
plane  wave.  If  the  incident  plane  wave  is  given  by 

^i (— t>  =  exp{-3k0i.xt}  (6-1-9) 

then  the  scattered  field  in  the  radiation  zone  of  the  tree  trunk 
has  the  form 

exp{- jk0x.  } 

ts(xt,u)  »  f(o,i)  ,  xt  =  jxj  (6-1-10) 

where  i  and  o  (=  x,t/xt )  are  unit  vectors  in  the  direction  of 
the  incident  wave  and  observation  point,  respectively.  Equation 
(6-1-10)  serves  as  the  defining  equation  for  the  scattering 
amplitude,  f(o,i_).  The  scattering  amplitude  is  directly  related 
to  the  Fourier  transform  of  the  transition  operator.  The  rela¬ 
tionship  has  been  derived  in  Appendix  B  of  [42]  and  is  given 
by 

f(o,i)  =  Yt(k0o,k0i,u>) 
where  t  is: 

t(kt,k£,u))  =  (2ir)  Jdx^dx^t  (xfc  ,x_ 

Although  the  transition  kernel  t(x,x',u>)  enters  into  the  defining 
equations  for  the  mean  intensity  and  the  correlation  function, 
it  is  the  scattering  amplitude  that  appears  in  the  final  results. 

Throughout  this  subsection,  the  w  dependence  of  quantities 
has  ,  been  explicitly  exhibited.  This  has  been  done  since  the 
calculation  of  the  space-frequency  correlation  function  to  follow 
requires  explicit  knowledge  of  the  frequency  dependence. 

6.1.2  Mean  Equation 

An  approximate  equation  for  the  mean  field  <i^>  was  obtained 
in  Section  2.2.3,  [42]  by  employing  the  Foldy-Lax  method  [19,35]. 
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The  method  is  valid  for  collections  of  two-dimensional  discrete 
scatterers  having  small  fractional  area.  This  criterion  is 
usually  satisfied  inasmuch  as  most  forests  have  fractional  areas 
ranging  from  .05  to  .2. 

The  approximate  equation  that  the  mean  field  satisfies 


(?£  +  k,)  <ip(xfc,u))  >  + Jd £tdx^.p(st)  t  (xfc  -  sfc ,x£  -  sfc , w)  <<Mx£,< 


=  -g(xt,ui) 


where  the  density  P(S£)  is  zero  for  s^.  outside  of  S,  i.e. 


(6-1-13) 


P(st)  = 


p(s.)  ,  s.  inside  S 


(6-1-14) 


,  sfc  outside  S 


This  equation  is  an  integro-dif f erential  equation  inside  the 
forest  but  when  s*.  is  outside  S,  the  mean  equation  reduces  to 
the  free  space  scalar  wave  equation. 

The  operator  notation  previously  introduced  can  be  used 
to  simplify  the  appearance  of  the  mean  equation.  By  using  the 
mean  operator  £  the  mean  equation  becomes 


£<ij/>  =  g 


(6-1-15) 


where 


£  =  L  -  y dstp (st)T(st) 


(6-1-16) 


Note  that  no  average  bar  appears  over  T  as  in  [42];  because 
for  cylindrical  scatterers,  T  *  T. 

The  mean  operator  £  has  been  introduced  since  it  is  important 
to  the  formulation  of  the  correlation  equation.  In  this  context 
its  inverse  will  often  be  required.  The  inverse  of  £  is  definted 
in  terms  of  the  mean  Green’s  function  G(xt,x£,w)  as 
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w)  <j>(x£,w) 


where 


£  G  =  6 


(6-1-17) 


(6-1-18) 


or  more  explicitly 


-<Vt 


+  k„)G(xt,xt 


/  n>) 


- Jd stdx^.p(st)  t(xfc 


st,oj)G(x^.,xt  ,w) 


= 6 (xt  -  xfc  )  (6-1-19 ) 

Thus  from  the  solution  of  Equation  (6-1-19)  the  inverse  of 
£“ 1  can  be  found  from  Equation  (6-1-17). 

In  the  low-frequency  limit,  the  transition  kernel  t  (xt,xt,u)) 
can  be  represented  as  Section  2.2.2,  [42] 

t(xt,x^,o))  =  k|a(a»)  6  (xt)  6  (x£)  (6-1-20) 

where  a  is  the  two-dimensional  polarizability  of  the  cylindrical 
scatterer.  For  a  circular  cylinder  of  susceptibility  X  it 
is  given  by 

* 

a(u))  =  -X^(u)27ra2  (6-1-21) 


where  a  is  the  radius  of  the  trunk.  If  the  low-frequency  form 
of  the  transition  kernel  as  it  appears  in  Equation  (6-1-20) 
is  used  in  Equation  (6-1-19),  the  equation  for  the  Green's  func¬ 
tion  simplifies  substantially.  It  is  given  by 

{ V2  +  k2([l  +  p  (xfc)  o  (w)3)G(xt,xt  ,w)  =  ~(5(£t-xt  )  (6-1-22) 

This  equation  will  be  employed  in  Section  6.3. 


6.1.3  Space-Frequency  Correlation  Equation 

The  basic  quantity  needed  to  characterize  communications, 
channel  is  the  space-frequency  correlation  function;  i.e.. 
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<ip(  5Ct ,  w)  ip*  (£t  •  >  •  This  represents  the  correlation  between 

a  field  component  of  frequency  w  at  position  x^  and  a  field 
component  of  frequency  ui  at  position  x^.  From  this  function, 
the  frequency  correlation  and  delay  spread  for  the  channel  can 
be  obtained  at  a  fixed  point  x^  ~  kt'  a°d  the  spatial  correlation 
can  be  found  at  a  fixed  frequency  u)  =  w. 

The  correlation  of  the  field  can  be  decomposed  into  a  co¬ 
herent  and  incoherent  part.  Consider  the  field  i(;(xt,u))  to  be 
represented  as  a  sum  of  its  mean  <ty(Xfbi)>  plus  its  fluctuating 
component  ij>(xt.,u))  so  that 

ip  ( x  t ,  a) )  =  <  tp  ( x  t ,  oj  )  >  +  i{;(xt,U))  (6-1-23) 

where  <  ipfxf- /x>)  >  =  0  by  definition.  By  using  this  equation,  the 
correlation  can  be  decomposed  as  follows 

<^(xt,u)iJ;*(xt/(D)>  =  - oj)  ><ip*  (xt,w)  > 

+  <iji(xt,w)Ij)*(xtrco)>  (6-1-24) 

The  first  term  <<|o<<p*>  is  the  mean  intensity.  It  represents 
the  complete  correlation  function  wheri  there  are  no  fluctuations. 
The  second  term  <<M*>  represents  the  effects  of  fluctuations. 
Thus,  the  representation  afforded  by  Equation  (6-1-23)  permits 
the  channel  to  be  divided  into  two  parallel  components  each 
with  its  own  distinctive  channel  characterization.  This  is 
shown  pictorially  in  Figure  6.2. 

The  equation  that  the  space-frequency  correlation  function 
satisfies  has  been  derived  in  Appendix  A.  The  methodology'  is 
the  same  as  the  derivation  appearing  in  Appendix  C  of  [42]  where 
the  space  correlation  equation  is  derived.  The  inclusion  of 
two  different  frequencies  only  changes  the  equation  slightly. 
The  equation  is  then  given  by 


where  tiie  notational  device  has  been  used  that  unstarred 


quantities  are 

quantities  are 

the  scattering 

functions  of  x*.  and  u>  while  starred  or  conjugated 
functions  of  xt  and  u .  The  integral  term  is 

or  diffusion  term  and  when  small 

££*<i(mJ>*>  5  gg* 

(6-1-26) 

or 

<M*>  =  £-1£*_1gg* 

(6-1-27) 

=  £~1g£*“1g* 

(6-1-28) 

■ 

(6-1-29) 

Thus  when 

the 

scattering  term  is  small. 

the 

random  component 

of  the  channel  can  be  neglected. 

In  the  following  subsections,  the  space-frequency  correlation 
equation  is  presented  €cr  the  two  special  cases  mentioned  earlier. 
First  in  Subsection  6.3,  the  frequency  correlation  function 
(:xt  -  xt)  is  given  for  an  infinite  forest  having  a  line  source 
excitation.  The  low-frequency  approximation  will  be  used  to 
make  the  equation  tractable.  In  Subsection  6.4  the  spatial 
correlation  function  (  w  *  w  )  is  given  for  a  half-space  forest 
with  a  normally  incident  plane  wave.  The  low-frequency 
approximation  is  not  employed  here,  so  the  results  are  valid 
for  all  frequencies  of  interest. 

6.2  Single-Scatterer  Characteristics 

The  properties  of  the  scatterer  ensembles  to  be  treated 
in  the  next  two  subsections  will  be  closely  related  to  the  char¬ 
acteristics  of  the  individual  scattering  elements.  In  this 
subsection,  the  differential  scattering  cross-section  and  total 
cross-section  of  an  isolated  tree  trunk  (circular  cylinder) 
are  plotted  for  typical  tree  trunk  diameters.  The  curves  obtained 
are  used  in.  latter  subsections  to  help  interpret  multiple  scat¬ 
tering  results. 


Consider  a  circular  lossy  dielectric  cylinder  of  radius 
a  whose  axis  coincides  with  the  z-axis.  The  cylinder  has  the 
complex  permittivity  -  jX£-  The  fundamental  quantity 

characterizing  the  scattering  properties  of  the  cylinder  is 
the  scattering  amplitude  as  defined  in  Equation  (6-1-10).  The 
scattering  amplitude  f  (o,  _i)  measures  the  response  of  the  scat- 
terer  in  the  direction  o  to  a  plane-wave  incident  upon  the 
scatterer  in  the  direction  i^.  These  unit  vectors  are  shown 
in  Figure  6.3  along  with  the  angle  of  incidence  4>^  and  the  scat¬ 
tering  angle  <$» s . 

The  dyadic  scattering  amplitude  for  a  cylinder  has  been 
computed  in  Section  3.1,  (42].  The  scalar  scattering  amplitude 

can  be  obtained  from  it  by  using  the  h°h°  component  of  the  dyadic 
at  Oi  =  tt/2.  The  result  is 


f(o,i)  =  fhh(0#i) 


0 . =n/2 

i 


V5 


j*/<  £  (.l)nchhei“«h-V 

n=-<® 


(6-2-1) 


The  expressions  for  the  Cn  are  given  in  Section  3.1,  [42]  and 
will  not  be  repeated  here.  The  are  functions  of  the  cylinder 
radius,  the  frequency,  and  the  complex  permittivity  X^. 

The  scattering  characteristics  of  the  cylinder  are  described 
by  the  differential  scattering  cross-section  od»  the  scattering 
cross-section  os,  the  total  cross-section  and  the  scattering 
albedo  w.  These  are  addressed  separately  in  the  subsections 
which  follow. 


6.2.1  Scattering  Formulas 

6 . 2 . 1 . 1  Differential  Scattering  Cross-Section 

The  differential  scattering  cross-section  is.  defined 


od(o,i)  »  |f(o,i)|  (6-2-2) 
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Using  Equation  (6-2-1)  in  the  above  equation  gives 


od(o,i)  =  od(4>*) 


too  -f-OO  £ 

-f-  E  E  (-D n+m  c*^  cjj11 

°  n=-°°  m=-» 


(6-2-3 


where  ip '  =  <Pi  -  4>s .  Thus  od  can  be  computed  from  a  knowledg 

of  the  Cn  and  it  is  only  a  function  of  the  difference  betwee 
the  incidence  angle  and  the  scattering  angle  $s.  The  differ 
ential  scattering  cross-section  also  has  the  following  tw 
properties : 


od(*')  =  odU*  +  2  it) 


(6-2-4 


od(y)  =  od(-$*)  (6-2-5 

Thus  it  is  seen  that  fljit')  is  an  even  periodic  function  o 

A  quantity  closely  related  to  o  d  is  the  bistatic  scatterin 
cross-section,  °bi»  which  is  given  by 

°bi  <°'±>  m  4T,°d<£»i.>  (6-2-6 

It  is  often  used  in  radar  applications.  The  differential  scat 
tering  cross-section  gives  the  gain  pattern  or  angular  distri 
bution  of  scattered  power  of  the  scatterer.  It  is  usually 
function  of  incident  angle  but  since  the  cylinder  is  rotationall 
invariant,  od  depends  only  on  the  difference  between  <j)  s  an 
♦i. 

6 . 2 . 1 . 2  Scattering  Cross-Section 

The  next  quantity  of  interest  is  the  total  power  scat 
tered  at  all  angles  surrounding  the  cylinder.  The  cross-sectio 
of  a  particle  which  would  produce  this  amount  of  scatterin 
is  called  the  scattering  cross-section  og  and  is  given  by 
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(6-2-7) 


/2ir  /.  2 TT 

°d(°-i>d*s  -J„ 


If  Equation  (6-2-3)  is  used  in  the  above  equation;  the  result 


4  v->  |  hi 

“s  =  T:  „?JCn 
n=-°° 


(6-2-8) 


6 . 2 . 1 . 3  Total  Cross-Section 

The  total  cross-section,  ot,  is  the  sum  of  the  scattering 
cross-section  os  and  the  absorption  cross-section  oa.  It  measures 
the  amount  of  power  the  scatterer  extracts  from  the  incident 
wave  in  scattering  and  absorption.  Since  oa  is  not  directly 
related  to  the  scattering  amplitude,  the  calculation  of  a ^  in 
terms  of  the  C^h  can  most  easily  be  done  by  using/optical  theorem. 
According  to  Appendix  F  of  142],  the  optical  theorem  for  two- 
dimensional  media  states 


o  =  -  Re[f  (i,i)e~jir/4] 

t  rr—  u -  J 


(6-2-9) 


Now  using  Equation  (6-2-1)  in  the  above  equation  yields 


4  ^  fhh"l 

°t  *  Rncn  J 


(6-2-10 ) 


6.2.1 .4  Albedo 

The  final  quantity  of  interest  is  the  albedo  W  which 
is  the  ratio  of  os  to  ot,  i.e. 


w  -  °s/ot 


( 6-2-11 ) 


The  albedo  always  lies  between  zero  and  unity  and  measures  the 
ratio  of  power  scattered  to  power  absorbed.  Scatterers  having 
albedos  close  to  zero  are  very  absorptive  while  scatterers  with 
albedos  close  to  one  scatter  much  more  power  than  they  absorb. 
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Single-Scatterer  Properties 


The  parameters  that  affect  the  cross-sections  are  the 
trunk  radius,  a,  the  frequency,  f,  and  the  complex 

susceptibility,  X^  .  The  scattering  cross-sections  will  be 

computed  for  trunk  radii  of  1,  5  and  10  centimeters  since  these 
represent  typical  trunk  radii  found  in  forests.  The  frequency 
range  will  be  chosen  from  200  MHz  to  2000  MHz  and  three  models 
for  the  susceptibility  will  be  used.  These  three  models  have 
been  called  CyberCom  Model  I,  II  and  III.  They  have  been  defined 
and  discussed  in  Appendix  A,  [42].  All  three  models  have  a 
constant  real  susceptibility  independent  of  frequency.  The 
imaginary  part  of  the  susceptibility  XT  differs  in  the  models. 

,  X# 

For  Model  I,  x£  is  constant;  for  Model  II,  the  conductivity 
of  the  medium  is  constant,  while  for  Model  III  relaxation  losses 
have  been  added  to  Model  II.  The  variation  of  X^  for  each  model 
is  shown  in  Figure  6.4.  Model  III  most  closely  depicts  the 
actual  trunk  medium.  The  ether  models  are  used  to  show  the 
sensitivity  of  the  results  to  variations  in  the  trunk 
susceptibility. 

6. 2. 2.1  Differential  Scattering  Cross-Section 

The  first  curves  to  be  discussed  are,  the  plots  of  o 
versus  f  shown  in  Figures  6.5,  6.6  and  6.7.  Each  figure  presents 
od  curves  (scattering  patterns)  at  frequencies  of  200,  300, 
and  600  MHz  at  a  fixed  tree  radius.  The  curves  have  been 
normalized  but  the  same  normalization  constant  is  used  for  all 
plots  having  the  same  tree  radius,  i.e.,  curves  on  the  same 
figure  have  the  same  normalization.  CyberCom  Model  III  has 
been  used  in  the  computations.  Scattering  patterns  have  not 
been  plotted  for  CyberCom  Model  I  and  II.  . 

An  examination  of  the  curves  shows  that  at  low  frequencies, 
i.e.  200  and  300  MHz,  the  trunks  of  radii  1  and  5  cm  are 
almost  perfect  isotropic  radiators.  As  expected,  as  the  frequency 
increases  for  a  fixed  trunk  radius,  the  scattering  pattern  becomes 
more  directed  toward  the  forward  direction.  This  trend  becomes 
more  pronounced  as  the  trunk  radius  is  increased w  . 


600  MHz 


Figure  6-5  Plane -Wave  Scattering  Pattern  From  a  Single  Trunk 
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Figure  6-6  Plane-Wave  Scattering  Pattern  From  a  Single  Trunk 
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6. 2. 2. 2  Albedo  and  Total  Cross-Section 

In  Figures  6.8  through  6.13  the  albedos  and  the  total 
cross-sections  are  plotted  as  functions  of  frequency  for  CyberCom 
Models  I,  II  and  III.  This  is  done  for  trunk  radii  of  1,  5 

and  10  centimeters.  The  albedo  curves  indicate  the  importance 

of  scattering  compared  to  absorption  while  the  a  t  curves  are 
directly  related  to  the  mean  wave  since  the  mean  wa''e  decays 
at  a  rate  exp[-patz).  Values  for  cg  have  not  been  plotted'  since 
they  can  be  computed  by  inrpection  from  the  albedo  and  curves 
by  using  Equation  (6-2-11). 

An  examination  of  the  albedo  graphs  shows  similar  curves 

for  all  permittivity  models  and  that  as  the  trunk  radius  becomes 

larger,  the  curves  oscillate  more  rapidly  with  smaller  amplitudes. 
This  is  because  large  trunks  reach  resonance  at  lower  frequencies. 
It  is  interesting  to  compare  the  ot  curves  for  1  and  5 

centimeters.  The  one  centimeter  curves  increase  with  increasing 
frequency  up  to  600  MHz.  This  indicates  the  attenuation  of 
the  mean  wave  is  becoming  greater  with  increasing  frequency. 
For  the  5  centimeter  trunks,  however,  the  o^.  curves  are  decreasing 
with  increasing  frequency  and  thus  the  mean  wave  propagates 
farther  as  the  frequency,  increases.  ,  , 

It  is  also  interesting  to  compare  Models  I,  II  and  III 
to  see  the  effects  of  loss  on  the  albedo  curves.  Examining 
Figure  6.10,  it  is  seen  at  the  high  frequency  part  of  the  band 
(1  GHz  <  f  <  2  GHz)  that  the  more  lossy  the  trunk  medium,  the 
smaller  are  the  oscillations  or  resonances  that,  occur.  It  should 
be  pointed  out  however,  that  these  oscillations  will  be  averaged 
out  when  a  distribution  of  trunk  sizes  is  used  to  represent 
a  realistic  forest. 

The  curves  of  the  albedo  for  trunk  radii  of  5  and  10  centi¬ 
meters  are  very  high,  indicating  that  the  scattered  component 
of  field  will  be  important.  Only  for  trunk  radii  of  1  centimeter 
and  for  frequencies  of  less  than  400  MHz  does  the  albedo  drop 
to  lower  values.  Thus  even  though  the  mean  wave  may  have 
decreasing  attenuation  for  increasing  frequencies  in  large  tree 
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Figure  6-10:  Albedo  Versus  Frequency  for  a  Trunk  Radius  of  5cro 
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radii,  the  scattered  component  may  be  much  larger  than  the  mean, 
making  the  behavior  of  the  mean  unimportant. 

6.3  Frequency  Correlation  Function  for  a  Thin-Trunk  Forest 


Propagation  of  wideband  signals  through  a  forested  environ¬ 
ment  requires  a  knowledge  of  the  frequency  correlation  function 
of  the  medium.  With  the  channel  characterized  in  this  manner, 
the  deleterious  effects  of  multiple  scattering  between  forest 
components  can  be  properly  assessed.  The  frequency  correlation 
function  can  be  found  by  solving  the  correlation  Equation  (6-1-25) 
with  jct  =  xt.  The  equation  is  difficult  tc  solve  unless  special 
assumptions  are  made.  In  the  present  case,  trunks  whose  radii, 
are  small  compared  to  a  wavelength  will  be  assumed. 

In  the  past,  several  investigators  have  found  approximate 
expressions  for  the  frequency  correlation  functioh.  This  calcu¬ 
lation  has  been  made  for  channels  having  fluctuating  media  such 
as  those  channels  involving  ionospheric  ( 41 , 53 , 77 , 78 ]  or  tropo¬ 
spheric  (34,35,691  turbulence,  as  well  as  channels  involving 
hydrometeors  C  3 5 1 .  In  all  of  these  applications,  the  forward 
scattering  assumption  has  ,  been  made  because  wavelengths  of 
interest  were  large  compared  to  medium  correlation  lengths  or 
particle  sizes.  In  the  present  case  of  forest  propagation, 
the  forward  scattering  assumption  is  not  valid  over  the  frequency 
band  of  interest  for  most  scatterers.  Thus,  the  correlation 
equation  must  be  solved  directly. 

In  this  subsection,  an  infinite  forest  is  considered  having 
constant  density  p.  The  waves  in  the  forest  are  excited  by 
an  infinite  line  source  located  at  the  orig.  - ;  thus  J*(xt)  * 
6(xt).  The  space-frequency  correlation  equation,  (6-1-25)  with 
xt  -  xt  will  be  solved  for  trunks  with  diameters  small  compared 
to  a  wavelength.  Using  this  low-frequency  assumption,  the 
correlation  equation  can  be  substantially  simplified.  To  see 
this,  multiply  Equation  (6-1-25)  from  the  left  by  and 
use  Equation  (6-1-20).  The  resulting  equation  for  the  frequency 
correlation  function  is; 
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i  (  Xt  t  i* '  i  t  ^  )  ~  X  ^  #  ^  1  9  ^  2  ^ 


-  z  }  Jc i£t  G  ( xfc  ~  £ fc ,  u  j )  G *  ( x t  -  sfc ,  w2 )  r  (sfc ,0)  J , ( 


(6-3-1) 


where  r  (it  •  ui  *<*;  *  =  < v (it  • w (xt  •  w2  ) >  is  the  frequency  cor¬ 
relation  function;  r  o  ( xt  ,cj  t  ,o>2  )  =  <  i|>(  xt  ,Wl  )  ><  tp*  ( xt  /(l>2  )  >  is  the 

coherent  portion  of  the  correlation  function. 


"  Bi?  =  pk\ ' 


(6-3-2) 


,k .  =  w./c. 

l  i  ' 


a .  =  a  ( ui .  ) 
1  x 


i=l  /  2 


(6-3-3) 


and  Co  -  (UoCo)_lj  is  the  free-space  velocity  of  light. 


The  Green’s  function  Gtx^w)  appearing  in  Equation  (6-3-1) 
satisfies  Equation  (6-1-9)  with  the  low— frequency  assumption 
and  *t0  ~  0.  The  equation  is 


(V2  +  k2 (1  + pa) )G(xt , w)  =  5 (xfc) 


(6-3-4) 


The  solution  to  this  two-dimensional  mean  Green's  function 
equation  is  given  by 


G(xt,u>)  «  j  Hi?)  (<xt) 


Xt  =  ‘^t1 


(6-3-5) 


where 


x  =  k/r  +  pn  ,  k  a  w/c, 


(6-3-6) 


and  h(*)(M  is  the  Hankel  function  of  zero  order  and  second 
kind.  The  mean  wave  <iHxt,u)>  can  easily  be  found  by  multiplying 
Equation  (6-3-4)  by  jwy0.  Referring  to  Equation  (6-1-13)  it 
is  seen  that 


<(Mxt,w)>  »  G(xt,w)/jwp0. 


(►  7) 


The  frequency  correlation  equation,  (6-3-1)  can  be  recast 
by  employing  the  Fourier  transform  since  the  integral  term  is 


4 
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a  convolution  of  GG*  and  f  .  Rather  than  proceeding  this  way, 
the  Green’s  functions  and  F0  will  be  replaced  by  their  far-field 
approximations.  This  simplifies  the  calculation  substantially 
and  only  restricts  the  observation  point  to  be  in  the  far-field 
of  the  source.  Tree  trunks  that  are  in  the  near  field  of  each 
other  are  not  correctly  accounted  for.  This  does  not  represent 
a  large  error,  however,  since  the  fractional  area  occupied  by 
tree  trunks  is  small. 

Now  using  the  asymptotic  approximation  for  the  Hankel 
function  in  Equation  (6-3-5),  the  Green's  function  becomes 


where 


G (x  ,  )  ~  J . .  e 

2/2nkxt 


-j(kxt  -kpy-  J) 


k  ~  k (1  +  p  a/2) 


(6-3-8) 


(6-3-9) 


has  been  used.  The  latter  is  a  good  approximation  for  k  since 
pa<<l  for  sparse  forests. 

If  this  far-field  approximation  is  used  in  the  frequency 
correlation  equation,  for  the  Green’s  function  and  for  r0»  the 
result  is  the  following  approximate  correlation  equation: 


r(xt,0Ji,w2)  =  K„  (xt,w,  ,u2) 

-  6, 2  J  dstK(xt  -  st,w-j  ,u>2)  T  ( x  t ,  tp  t  ,u>2) 


(6-3-10) 


where! 


with 


c  ”j7Xt 

K(xt,u),,o,2)  =  —  e 


(6-3-11) 


-  Ak  -  ptkjO,  -k2a*)/2  ,  Ak  =  kj-k. 

c  —  ((8n)  2k1k2rl* 


(6-3-12) 


(6-3-13) 
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» 71  n  v-  ^  *c-»  o»  •  « •  /  * 


K0  (xfc  ,u1  ,m2)  =  K  ( x  t ,  oj  x  ,(02)/(u>mo) 


(6-3-14 


This  approximate  frequency  correlation  equation  is  easier  t 
transform  since  the  Fourier  transform  of  K(x^,^ltu)2)  can  b 
calculated  explicitly. 


Taking  the  Fourier  transform 


f(k.,ui'.,td ,)  =  — - —  fdx J(x  ,m  ,io  ) 
~t  (2  _t  -*■ 


(6-3-15 


of  the  correlation  Equation  (6-3-10)  and  solving  for  f (k^ fWj ,w2 
yields 


n  k  t ' w  i ' 00  2 )  = 


Ko 

l-6i2K(kt,avu>2) 


,  B12=(2r)2B1 


(6-3-16 


where 


K  ( k  t ,  w  j  ,  a)  2 )  = 


2ttA£  +yz 


kt  »  |ktl 


'o  ( k t , a> J  ,W2)  =  K(kt,Wj  ,u)2)/tu2y; 


(6-3-17 


(6-3-18 


Y  =  3Y 


(6-3-19 


The  solution  for  the  frequency  correlation  function  is  not 
obtained  by  transforming  Equation  (6-3-16).  This  calculatioi 
yields 


/•«*  dk  k.  J0  (k  x  ) 

:x.,u)  .u  )  -  cl.  — - 

1  1  2  V  a?  ;  ?2  - 1 


(6-3-20 : 


where 


H..,c  3/ 

b  =  ■  -  =  p(k,kj  /2  a, a*/ 4 


(6-3-211 


In  obtaining  Equation  (6-3-20)  the  integration  over  the  angular 
variable  <t>  has  been  performed  giving  the  Bessel  function  J0(ktxt)* 
The  calculations  of  this  subsection  are  similar  to  those  appearing 
in  Chapter  7,  [42];  in  fact,  this  calculation  reduces  to  that 

one  for  uj  =  io2  =  to. 

The  integral  representation  for  r(xt,(jjj  ,to2 )  is  very  slowly 
convergent  and  thus  difficult  to  numerically  evaluate.  Following 
a  parallel  development  to  that  performed  in  Appendix  G,  [42], 
an  alternate  representation  for  F  ,oo2  )  can  be  obtained. 

It  is 

F  (xfc,to1  ,to2)  =  -  jncbH  j  ^  (-j/y2  “  *>2  xfc) 

f -  H<2)(- 

-  j  cl  dp  p2 - - - —  (6-3-22) 

•^0  p2  +  b2 

Before  proceeding  to  evaluate  numerically  the  above 
expression  it  is  useful  to  represent  the  correlation  function 
as  a  sum  of  a  coherent  part  T0  and  an  incoherent  part  Fj,  i.e. 

r(xt,u)1,u)'2)  =  F0{xt,wl  ,u>2)  +  Fx  (xt,toi  ,to2)  (6-3-23) 

It  is  the  ihcoherent  portion  of  the  correlation  function  Tj 
that  should  tend  to  zero  as  |u>j-u)2|  becomes  large. 

The  incoherent  correlation  Fi  can  be  obtained  by  subtracting 
To  from  r.  An  expression  for  To  can  be  obtained  from  Equation 
(6-3-22)  by  setting  S12  =  0,  which  implies  that  b  *  0.  The 
result  is 

00  '  _ ' _  1 

ro(xt,wi,u)2)  =  -jcj  dp  rti  2  ^  (-j/p2  +  y 2  xt)  (6-3-24) 

.-'o’ 

Now  subtracting  To  from  T  as  given  in  Equation  (6-3-22)  yields 
the  desired  expression  for  Tj 


I’j  (xt,(i)1  ,U2)  =  -  j iTcbHo  (-j/y2  -  b2  xt) 


rm  (-j/p2  +  y  x.  ) 

-  jc  I  dp  p2 -  (6-3-25 

J0  p2  +b2 

A  normalized  form  of  Tj  has  been  plotted  in  Figures  6.1 
and  6.15  for  center  frequencies  of  200  and  300  MHz  respectively 
For  each  frequency  a  plot  is  shown  for  distances  of  200  ant 
500  meters.  Only  frequencies  close  to  the  lower  edge  of  th< 
band  of  interest  and  trunks  of  small  radius  are  considered  becaust 
of  the  Rayleigh  assumption  implicit  in  the  derivation  of  Equatioi 
(6-3-1).  The  normalized  correlation  function  Tj  represent: 
divided  by  its  value  at  Wj  =  iu2,  i.e. 

C j  (x^ / w j  ^ u>2 )  =  r  j  ( i w j  ,w2)/ri( xfc , uj j  , u) j )  (6—3  —  26 

Thj.s  has  been  plotted  versus  the  frequency  difference  -w2 
Both  the  real  and  imaginary  parts  of  have  been  plotted  ,  ii 

each  case.  The  upper  and  lower  envelopes  represent  graphs  ol 
the  magnitude  of  Tj  and  its  negative  respectively.  The  curves 
are  plotted  until  the  value  of  | T^x  I  reaches  0.1  and  then  the} 
are  terminated.  The  correlation  bandwidth  Bc  measures  tht 
frequency  difference,  w,  -  w2  for,  which  Tj  is  negligible.  II 
will  be  assumed  that  Bc  is  the  value  at  which  j  T  j  |  is  equaJ 
to  0.1. 

The  curves  will  now  be  discussed  wholly  in  terms  of  tht 
correlation  bandwidth,  Bc.  Examining  Figure  6.14  for  200  MHz, 
it  is  seen  that  the  correlation  bandwidth  decreases  from  6  MHa 
to  4  MHz  as  the  distance  increases  from  200  to  500  meters. 
Increasing  the  frequency  by  50%  to  300  MHz  lowers  the  correlatior 
bandwidth  to  2  MHz  at  200  meters  and  1  MHz  at  500  meters.  The 
trends  appear  to  be  that  the  correlation  bandwidth  decreases 
with  increasing  distance  and  increasing  frequency.  It  snoulc 

be  recalled,  however,  that  this  is  just  a  low-frequency  result, 

> 

and  the  trends  should  not  be  interpreted  too  broadly. 


y  -  200.  MHz 

no  =  1000. /ha 


Rad  -  * 


D  i  st  an 


Distance  «  500 

Tau  "  .83 


Versus^DifferenceeprequencynCt^°nS 


for  a  Forest  of  Resonant  Tree  Trunks 


Propagation  through  a  random  medium  of  discrete  scatterers 
such  as  a  forest  is  characterized  by  a  two-component  channel 
model.  This  model  is  illustrated  in  Figure  6.2  and  consists 
of  mean  and  random  or  fluctuating  components.  In  this  subsection, 
the.  properties  of  these  components  will  be  studied  in  a  trunk- 
dominated  forest  of  resonant  trees.  Emphasis  will  be  placed 
on  obtaining  the  two-point  correlation  function  <ip  (x)  ip*  (jc)  > 
at  one  frequency  u> .  Two  properties  of  major  interest  are:  the 
relationship  of  coherent  to  incoherent  intensity  fluctuations 
as  a  function  of  distance  and  frequency;  and  the  behavior  of 
the  space  correlation  function,  <ip  (x )  <p  (x)  > ,  as  the  distance 
between  x  and  x  increases. 

The  work  of  this  subsection  in  part  represents  a  continuation 
of  the  material  in  Chapter  7,  [42].  There  the  relationship 

between  the  coherent  and  incoherent  intensity  has  been  examined 
for  an  infinite  forest  of  thin  tree  trunks.  The  thin  tree  trunk 
assumption,  although  valid  for  only  the  lower  portion  of  the 
frequency  band  of  interest,  allowed  the  correlation  Equation 
(6-1-25)  with  ?!  -  x  and  u  =  u  to  be  solved  exactly.  The  results 
showed  that  for  frequencies  as  low  as  300  MHz  intensity  fluctua¬ 
tions  are  quite  important  and  that  they  increase  with  distance. 
These  findings  will  now  be  extended  to  the  case  of  resonant 
trunks. 

For  the  resonant  trunk  case,  solution  of  the  correlation 
equatidn  by  a  transform  technique  does  not  seem  possible,  instead, 
a  two-variable  perturbation  procedure  is  used  to  derive  a  gener¬ 
alized  transport  equation  from  the  correlation  equation.  The 
fractional  area  is  the  small  perturbation  parameter  used  in 
the  above  analysis.  The  derivation  of  the  transport  equation 
from  the  correlation  equation  is  presented  in  Appendix  C.  The 
analysis  of  the  transport  equation  will  be  given  in  this 
subsection.  Transport  theory  has  been  applied  extensively  to 


discrete  scatter  problems  [25,34,54,58],  however,  CyberCom 

has  found  nc  application  to  tree  trunk-dominated  forests. 

6.4.1  Problem  Formulation 

Consider  a  trunk-dominated  forest  that  occupies  the  half 
space,  z  >  0,  as  is  shown  in  Figure  6.16.  The  forest  is  composed 
of  parallel  circular  dielectric  cylinders  having  radius  a.  They 
are.  distributed  with  constant  density  p.  A  plane  wave  of  unit 
intensity  is  normally  incident  on  the  half  space  and  the  space 
correlation  function  <  ip(  x )  i|>(>c)  >  is  evaluated  at  x  and  x  as  shown 
in  Figure  6-16.  To  allow  for  certain  simplifications  due  to 
symmetry,  the  two  points  will  be  assumed  to  be  in  the  same 
transverse  plane,  i.e.,  their  z  coordinates  will  be  the  same. 
This  restriction  will  limit  the  investigation  to  transverse 
correlations,  however,  this  is  the  case  of  most  practical  interest 
for  diversity  antenna  siting. 

Because  of  the  assumed  location  of  the  observation  points 
x  and  x  and  the  uniformity  of  the  problem  in  the  transverse 
direction,  the  correlation  function  T(x,x)  =  <i|>(x)  (x)  >  only 

depends  on  the  difference  of  the  transverse  coordinates  x-x 
and  the  normal  distance  into  the  half  space  z,  thus 


r  =  nz,rx) 


rx  =  X  -  X 


(6-4-1) 


This  correlation  function  can  be  written  in  terms  of  a 
generalized  specific  intensity  function,  J(z,rx, $)/  as  follows: 


r(2-rx> 


ik0r  sin$ 

d<j>  J  (z,rx,$)e  x 


(6-4-2) 


When  rx  =  0  the  generalized  specific  intensity  reduces  to  the 
specific  intensity  I(z,<p),  i.e. 


J  (z,0,  <p)  =  I(z,<f>) 


(6-4-3) 


From  Appendix  C,  it  is  seen  that  J(z,rx,  )  satisfies  the  following 
transport  equation: 


cos  <J> 


dJ(2,rx,4>) 


d T - +PatJ<z,rx,c|>) 

(6-4-4  ) 

/•2ir 

xsin<J>)l  d<{),obi(4),d)' )  I  (z,^' ) 

<P  e  [0 , 2 rr] 


=  pA(r 


where 


A(£)  =  exp{i2n2pf (o,o) ^/yk0 } 


(6-4-5) 


,-v 

■> 

iU 

« 


Here  A ( c )  is  the  slowly  varying  part  of  the  mean  Green's  function 
discussed  in  Appendix  B.  The  forward  scattering  amplitude  f (o,o) 
appears  in  Equation  (6-4-5).  It  is  in  the  direction  jo,  a  unit 
vector  making  an  angle  <j>  with  the  z-axis  as  in  Figure  6-17. 

The  solution  of  Equation  (6-4-4)  for  the  generalized  specific 
intensity  when  used  in  Equation  (6-4-2)  will  give  the  desired 
space  correlation  function.  This  equation  will  be  solved  in 
two  stages.  First,  rx  will  be  set  equal  to  zero  and  the  equation 
will  be  solved  for  the  specific  intensity  I(z ,<}>);  second,  the 
solution  for  J(z,rx,<J>)  will  be  found  in  terms  of  I(z,<J>).  It 
should  be  pointed  out  that  the  transport  equation  for  J(z,rx,<J>) 
is  the  standard  radiative  transfer  equation  that  has  been 
considered  by  many  other  research  workers  when  rx  is  set  equal 
to  zero. 

6.4.2  Intensity  Fluctuations 

In  this  subsection,  the  transport  equation  for  the  specific 
intensity  will  be  solved  numerically  and  the  results  discussed. 
The  transport  equation  for  the  specific  intensity  is  given  by 


COS  (j)  : 


:dl(z,4>) 


dz 


/•2tt 

+  potI(z,<f>)  -  pj^  obiU  -  $')  I(z,V)d*' 


(6-4-6) 


where  the  rotational  symmetry  of  the  cylindrical  scatterers 
has  been  utilized  by  replacing  a^i (  <f>,  <j>'  )  by  Obi (<}> - <J>'  )  (see 
Equation  (6-2-3).  To  solve  the  transport  equation,  the  boundary 
conditions  at  z  «  0  and  z  *  ®  will  be  needed.  These  boundary 
conditions,  are: 
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*  •  V  V  1 


1(0, (j))  =  6(4)) 


I 


<p  £  <t> 


(6-4-7) 


=  o  ,  <J>  e  (0,271]  (6-4-8) 

where  $+  is  the  set  of  angles  <|>  corresponding  to  propagation 
in  the  forward  direction.  The  first  boundary  condition  as  given 
by  Equation  (6-4-7)  states  that  the  incident  energy  on  the  half 
space  is  in  the  <£  =  0  direction  only.  The  second  condition 
says  that  all  energy  will  be  absorbed  as  z  approaches  infinity. 

The  specific  intensity  is  now  decomposed  into  two  parts 
called  the  reduced  specific  intensity  Ir(z,<J>),  and  the  diffuse 
specific  intensity  Id(z,$).  Thus 

I(z,4>)  =  Ir(z,<}>)  +  Id(z,<J>)  (6-4-9) 

where  Ir  is  the  solution  of  the  transport  equation  without  the 
scattering  term.  It  satisfies  the  equation 

dl  (z,$j 

cos  <J> - ^ -  +  potIr(z,c|>)  =  0  (6-4-10) 

with  boundary  conditions  as  given  in  Equations  (fc-4-7)  and 
(6-4-8).  The  equation  for  the  diffuse  intensity  is  obtaineo 
by  subtracting  Equation  (6-4-10)  from  Equation  .(6-4-6).  The 
result  is 

dld  (z ,  <p)  /•  2tt 

cos  (J> - — +  potId(z,4>)  =  obiU  -  <J>')  Id(z,<j>,)d<J>* 

Z-2  TT 

+  pjo  °bi(^  -  <*>* )  Ir(z,4>»  )d(J>»  (6-4-11) 

with  boundary  conditions 

Id«M)  =»  0  ,  (j)  e  $+  (6-4-12) 


*  0 


9 


$  e  [0,2ir] 


(6-4-13) 
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It  should  be  pointed  out  that  the  reduced  and  diffuse  specific 
intensities  are  exactly  the  same  as  the  mean  or  coherent  and 
the  incoherent  specific  intensities  respectively.  These  quanti¬ 
ties  when  substituted  in  Equation  (6-4-2)  yield  the  coherent 
and  incoherent  intensities. 


The  solution  of  the  reduced  equation  is  obtained  by  inspec¬ 
tion.  It  is 


I  (z,<fr)  =  6(4>)e 


-T 


,  t  =  patz 


(6-4-14) 


where  t  is  the  optical  depth.  Next,  the  reduced  intensity  is 
substituted  into  the  right-hand  side  of  the  diffuse  Equation 
(6-4-11).  Before  solving  this  equation,  the  symmetry  of  I(z,<J>) 
about  <{>  =  7T,  i.e., 


I  (z  ,  <j>)  =  I  (z ,  2ir  -  <{>) 


<t>  e  E  0  ,  tt  ) 


(6-4-15) 


is  used  to  reduce  Equation  (6-4-11)  from  the  angular  range 
cj>e  [  0 , 2  tt]  to  the  angular  range  <f>e[0,ir].  The  diffuse  equation 
becomes 


and 


dld(x,y) 


■  + 


Id(T,p)  =  f  *  Id(x,yk)dy’ 

•'-l  /I  -  y  ’  2 


(6-4-16) 


+  !bi^ie-x 


ye  [-1,1] 


ld(0,y)  =  0 


y  e  [0,1] 


(6-4-17) 


id(°°/  y)  =  o 


y  e  [-1,+1] 


(6-4-18) 


The  above  equation  has  been  written  in  terms  of  the  dependent 
variable  y  (y  =  cos  <j>  )  and  x  instead  of  $  and  z.  The  phase 
function  p(y,'i')  is  given  by 


c,.  (<j>  +  $ ' )  +  a.  .($-$') 

p  ( y ,  y ' )  =  — - - - — - 

°t 


(6-4-19) 


The  integro-dif ferential  equation  fqr  l£<x,y)  cannot  be 
solved  analytically  unless  o  bi  is  independent  of  angle  (isotropic 


6-39 


scatterer)  [20]  which  is  not  the  present  case  as  evidenced  by 
the  radiation  patterns  shown  in  Figures  6.5  thru  6.7.  Numerical 
methods  must  be  used  to  obtain  quantitative  results.  The  method 
of  discrete  ordinates  [20]  will  be  employed.  The  integral  term 
is  represented  by  a  discrete  sum  using  the  Chebyshev-Gauss  quad¬ 
rature  formula 


J.  /I  -  n=-N  J 

_1  j^O 


(6-4-20) 


where  the  Mj's  correspond  to  discrete  angles  <j>j  (<f>j  =  cos-1yj) 
which  are  defined  by 

!sin  (2  j  -  1)  tt/4N  j=l,...,N 

(6-4-21) 

M  ,  j— 1,...,-N 


Using  the  quadrature  formula  in  the  transport  equation, 
(6-4-16)  to  approximate  the  integral,  and  evaluating  the  equation 
at  y  =  Pj,  j  =  _+l,...,_+N,  results  in  a  system  of  2N  ordinary 
differential  equations  with  constant  coefficients.  These 


equations  are: 


+  T77  d  (T^Mi}  =  2N  V T 


Id(T,Uj) 


+  -Hi - i_  e  T  ,  i=+l,...,+N  (6-4-22) 

°t  ~ 


with  boundary  conditions 


ld(0.,ui)  =  0 


i=l, . . . ,N 


(6-4-23) 


Id^a>,yi)  =  0  *  i=+l,...,+N 


(6-4-24) 


The  equations  can  be  solved  for  the  Id(  ) ,  i  =  +1,...,+N, 
by  standard  numerical  techniques.  The  larger  the  value  of  N 
chosen,  the  closer  the  computed  values  come  to  the  actual  solution 


of  Equation  (6-4-22),  and  the  more  points  are  available'  to  approx¬ 
imate  1,3  (  t,u)  over  the  complete  angular  range. 

The  calculation  of  I^CtiH)  has  been  carried  out  for  N  =  20. 
Ail  calculations  in  this  subsection  have  been  done  using  CyberCom 
Model  III  for  the  trunk  permittivity,  a  tree  radius  of  5  centi¬ 
meters  and  a  density  of  1000  trees  per  hectare.  Curves  of  the 
normalized  diffuse  specific  intensity,  1,3,  versus  the  azimuthal 
scattering  angle,  <j>,  have  been  plotted  for  frequencies  of  200, 
300,  and  600  MHz  and  optical  depths  of  0.5,  1  and  5.  This  is 

a  total  of  nine  graphs  which  have  been  included  in  the  Annex 
with  figure  numbers  6A-1  through  6A-9.  In  addition,  three  curves 
have  been  plotted  at  a  constant  distance  of  200  meters  for  fre¬ 
quencies  of  200,  300,  and  600  MHz.  These  are  included  in  the 
Annex  with  figure  numbers  6A-10  through  6A-12.  The  curves  have 
been  presented  in  a  polar  format.  A  spline  routine  has  been 
used  to  connect  the  discrete  output  results  of  the  computation 
by  a  smooth  curve. 

An  examination  of  the  complete  set  of  curves  shows  several 
trends.  The  first  trend  is  an  increase  in  the  amount  of.  forward 
scattering  as  the  optical  depth  increases  for  fixed  frequency. 
.This  effect  is  illustrated  in  Figures  6-18  and  6-19.  The  first 
figure  is  for  an  optical  depth  of  0.5  and  it  shows  a  relatively 
large,  amount  of  backscatter.  The  second  figure  is  for  an  optical 
depth  of  5  and  it  shows  that  the  radiation  is  almost  ,  all  in 
the  forward  direction. 

This  effect  has  been  observed  by  Whitman  et  al.  [75]  in 

their  examination  of  .millimeter  wave  scattering  from  a  half 

space  of  three-dimensional  scalar  isotropic  spheres.  To  aid 
in  making  a  comparison  with  their  work.  Figures  6A-1,  6A-2  and 
6A-3  have  been  redrawn  in  rectangular  format  and  presented  in 

Figure  6-20.  A  comparison  of  this  figure  with  several  of 

Schwering  et  al.  [58]  shows  quite  good  qualitative  agreement. 

The  second  trend  can  be  observed  by  examining  Figures  6A-10 
through  6A-12.  There,  for  a  constant  distance  'of  200  meters, 

the  amount  of  forward  scattering  increases  as  the  frequency 
increases.  This  increase  ir.  forward  scattering  can  be  explained 
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Figure  6-19i  Normalized  Specific  Intensity  Versus  Azimuthal  Scatter  Angl 


Figure  6-20 (b):  Normalized  Specific  Intensity  Versus  Azimuthal  Scattering  Angle 
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by  examining  the  radiation  pattern  of  the  scatterer  (Figure 
6.6).  There  it  is  observed  that  the  radiation  pattern  becomes 
more  forward-directed  as  frequency  increases. 

The  specific  intensity  at  a  point  represents  the  amount 
of  energy  passing  through  that  point  at  a  specific  angle.  When 
this  specific  intensity  is  summed  over  all  angles,  the  total 
intensity  at  a  point  is  obtained.  To  see  this,  rewrite  Equation 
(6-4-2)  for  rx  =  0.  The  result  is 

2  /*2u 

<U(X)|  >  =  I  ( z)  =  /  a#  I  (z,4>)  (6-4-25) 

J0 


where  the  correlation  function  r(z,0)  has  been  set  equal  to 
I(z),  i.e.,  r(z,0)  =  I(z). 


This  intensity  can  be  broken  up  into  a  mean  or  coherent 
component  IQ(z)  and  a  random  or  incoherent  component  Ij(z). 
These  components  correspond  to  angular  integrals  over  the  reduced 
and  diffuse  specific  intensities  respectively.  From  Equations 
( 6  —  4  —  25 )  and  (6-4-9),  it  is  found  that 


I0(z) 


(6-4-26) 


and 


IjU) 


(6-4-27) 


The  coherent  intensity  can  be  obtained  directly  by  using  result 
(6-4-14)  in  (6-4-26);  thus 

I0(z)  -  e"T  ,  t  =  potz  (6-4-28) 


The  angular  symmetry  of  the  diffuse  specific  intensity  about  <p  =  * 
can  be  used  to  reduce  the  integral  given  in  Equation  ( 6  —  4  —  27 )  to 
the  interval  [  0 , tt  ) ;  then,  a  transform  to  p  =  cos  <p,  results 
in  the  following  expression  for  Ij(z): 
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0 


(6-4-29) 


r +1  Ij(t,p) 

IT  (z)  =  2  I  — . .  dy 

-l  /I  ~  li2 

Now  employing  the  Chebyshev-Gauss  quadrature  formula  given 
in  Equation  (6-4-20),  an  approximate  expression  for  the  incoherent 
intensity  is  found  in  terms  of  the  numerical  results  previously 
obtained.  The  incoherent  intensity  is  expressed  as 

N 

II(z)  =  £  £  Id{T'Mi)  (6-4-30), 

j=-N  J 

j^O 

The  results  of  the  computer  calculations  will  be  plotted 
in  three  formats.  These  formats  are:  normalized  intensity 

versus  distance;  incoherent-to-coherent  intensity  ratio  versus 
distance;  and  finally,  normalized  coherent  and  normalized 
incoherent  versus  distance.  The  intensity  has  been  normalized 
to  the  Value  of  the  incident  intensity.  Curves  will  be  drawn 
for  frequencies  of  200,  300,  and  600  MHz  and  radii  of  1,  5, 

and  10  centimeters.  All  curves  have  been  included  in  the  Annex 

under  figure  numbers  ,  6A-13  through  6A-21.  Sample  curves  have 
been  placed  in  the  text  bearing  figure  numbers  6.21  and  6.22 
for  illustrative  purposes. 

An  examination  of  the  normalized  coherent  intensity  as 
shown  in  Figures  6A-15,  6A-18  and  6A-21  shows,  that  all  curves 
decrease  with  increasing  distances.  The  rate  of  decrease  is 
directly  related  to  a*,  since  the  normalized  coherent  intensity 
varies  as  exp(-potz).  This  can  be  verified  by  referring  to 
the  total  scattering  cross-section  curves  shewn  in  Figures  6.9, 

6.11  and  6.13.  (Note  tick  marks  have  been  included  at  300  and 
600  MHz.)  It  can  be  seen  that  both  the  coherent  intensity  values 
and,  the  ot  values  increase  with  frequency  for  a  trunk  radius 
of  1  centimeter,  while  they  both  decrease  with  frequency  for 
trunk  radii  of  5  and  10  centimeters. 

In  Figures  6A-14,  6A-17  and  6A-20  the  ratios  of  the  incoherent 
to  coherent  intensities  have  been  plotted.  The  basic  trend 
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these  curves  show  is  that  the  incoherent  wave  becomes  increasingly 
important  as  frequency  and  distance  increase.  Since  the  amount 
of  scattering  directly  affects  the  importance  of  the  diffuse 
term,  it  is  interesting  to  compare  the.  albedo  curves  presented 
in  Figures  6.8,  6.10  and  6.12,  to  the  curves  for  the  intensity 

ratio  under  consideration.  The  comparison  shows  that  in  almost 
all  cases  an  increase  or  decrease  of  albedo  corresponds  tp  an 
increase  or  decrease  in  the  incoherent -to  -  coherent  intensity 
ratio.  This  is  not  the  case  for  curves  having  a  trunk  radius 
of  10  centimeters  where  the  200  and  600  MHz  curves  are  inter¬ 
changed.  It  should  be  noted,  however,  that  the  albedo  values 
and  the  curves  are  close  to  each  other  in  this  case. 

6.4.3  Space  Correlation  Function 

The  space  correlation  function  T  (z,rx)  can  be  found  by 
solving  the  generalized  transport  Equation  (6-4-4)  for  J(z,rx,<j>) 
and  then  using  this  quantity  in  Equation  (6-4-2)  for  F(z,rx). 
Because  of  the  transverse  placement  of  the  observation  points 
x  and  x,  the,  generalized  specific  intensity  J(z,rx,<f>)  can  be 
found  directly  in  terms  of  I(z,<J>).  It  can  be  verified  by  direct 
substitution  that 


J(z,rx,<|>)  =  A(rxsin9)I(z,4)  (6-4-31) 


and  thus  J  can  be  obtained  by  using  the  intensity  values  already 
computed . 

Substituting  Equation  { 6—4—31 )  into  the  expression  for 
the  intensity  gives 


/•2tt 

»r„)  =  /  dd>A (r  s 


T  (z 


'I 


0 


d<J>A(r  sin6)I(z,<J>)e 


ik„r  sinp 


2 tt  i<r  sin<{! 


d4>I(z,<j>)e 


where 


k  =  k.  +  <5k 


(6-4-32) 

(6-4-33) 


(6-4-34) 
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I 


(6-4-35) 


6  k  =  6  k 


j  6  k  11  =  27T2pf  (o,o)/yk0 


Here  k  is  the  mean  wave  propagation  constant  given  in  Equation 
(4-2-30)  of  [42]  with  0^  =  0. 

Before  proceeding  further,  the  correlation  function  is 
broken  up  into  coherent  and  incoherent  components  as  follows: 

r(z,rx)  =.r0(z,rx)  +  ri(z,rx)  (6-4-36) 

where 

-P 

,  r0(z,rx)  =  e  r  (6-4-37) 

and 

^2 7r  ixr  sin4> 

^(z,^)  =  d<pid(z,<t>)e  X  (6-4-38) 

Here  r0(z,rx)  is  the  coherent  portion  of  the  correlation 
function.  It  should  be  noted  that  it  does  not  depend  on  rx; 
thus  the  mean  fields  at  x  and  x  remain  completely  correlated 
for  all  values  of  rx.  This  is  not  surprising  since  the  mean 
field  is  a  deterministic  quantity. 

\ 

The  incoherent  portion  of  the  correlation  function  is 
F  i  ( z , rx ) ,  Using  the  symmetry  of  Id(z,4>)  about  <{>  =  it  ,  Equation 
(6-4-38)  can  be  rewritten  as 


fir  -6ic"r  sin<|)  , 

Tj  ( z , rx)  -  2j  d4>Id(z,4i)ccs[  (k0  +  6K')rxsin$]e 

(6-4-39) 

It  can  be  seen  from  Equation  (6-4-39)  that  the  correlation 
function,  Tj,  i  r?  real.  This  is  not  a  general  property  of  Tj 
but  results  from  the  transverse  location  of  the  observation 
points  x  and  x. 

The  numerical  analysis  of  the  integral  for  Fj  must  be  treated 
with  care.  The  function  !<} ( z,  $)  and  exp(-<5<”rxsin  $  )  are  slowly 


varying  functions  of  angle,  while  the  cosine  function  becomes 
a  rapidly  varying  function  of  angle  as  k0rx  becomes  large. 

To  see  the  behavior  for  large  k0rx,  the  integral  can  be 
asymptotically  evaluated  by  the  method  of  stationary  phase. 
The  stationary  point  occurs  at  <p .  =  tt/2  and  the  resulting  asymp¬ 
totic  evaluation  yields 

rrr  kx| 

ri(z,rx)  2  yk  ■—  ■  Id(z,ir/2)e 

cos[(k0  +  5^")rx  +  7r/2]  (6-4-40) 

This  approximate  expression  shows  that  for  large  k0rx  the 
correlation  function  falls  off  at  the  decay  rate  of  the  mean 
wave.  The  validity  of  the  above  asymptotic  expansion  requires 
the  ,  I<j(z  ,<J> )  be  a  slowly  varying  function  of  <P .  As  z  becomes 
large,  I^i  ( z ,  4> )  varies  more  rapidly  with  <j> .  Thus  as  z  becomes 
large,  the  correctness  of  the  above  result  falls  into  question. 

The  correlation  function  has  been  plotted  for  frequency 
values  of  200,  300  and  600  MHz  in  Figures  6.23,  6.24  and  6.25 
for  distances  of  100  and  500  meters.  A  spline  curve  fit  to 
I<j(z,<J>)  has  been  used  to  numerically  evaluate  the  4>  integral. 
The  computed  correlation  function  in  the  figures  is  shown  by 
the  thin  solid  line.  The  thick  solid  line  is  the  envelope, 
and  the  dashed  line  is  the  asymptotic  expression  given  in  Equation 
(6-4-40). 

An  -  examination  of  the  correlation  curve  shows  that  for 
200  and  300  MHz  the  correlation  function  does  not  change  sub¬ 
stantially  as  the  distance  varies  from  100  to  500  meters.  It 
should  also  be  noted  that  at  these  frequencies,  the  asymptotic 
expansion  is  quite  close  to  the  cpmputed  curve  for  most  of  the 
plotted  range.  Referring  to  the  approximate  correlation  func¬ 
tion  given  in  Equation  (6-4-40),  it  is  seen  that  the  correlation 
distance  does  not  depend  on  rx.  This  explains  the  invariance 
the  correlation  length  with  distance  observed  in  Figures 
6.23  and  6.24. 
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The  situation  is  different  in  Figure  6.25  at  600  MHz.  Here 
the  la (z,$)  becomes  more  rapidly  varying  for  large  z  and  the 
asymptotic  expansion  is  not  valid  for  a  large  portion  of  th>- 
z  =  500  meter  curve.  Thus  it  is  seen  for  this  case  that  the 
space  correlation  function  appears  to  decrease  with  distance. 

As  a  final  observation,  it  should  be  noted  that  the  space 
correlation  length  decreases  with  increasing  frequency.  This 
is  the  expected  result. 
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Figure  6A-ls  Normalized  Specific  Intensity  Versus  Azimuthal  Scatter  Angles 
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Figure  6A-2:  Normalized  Specific  Intensity  Versus  Azimuthal  Scatter  Angles 
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Figure  6A*-8;  Normalized  Specific  Intensity  Versus  Azimuthal  Scatter  Angles 
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Figure  6A-11;  Normalized  Specific  Intensity  Versus  Azimuthal  Scatter  Angles 
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APPENDIX  A 


Space-Frequency  Correlation  Equation 

In  this  Appendix  the  approximate  equation  for  the  space- 
frequency  correlation  function  will  be  derived.  The  methodology 
employed  in  the  derivation  will  parallel  the  Foldy  development 
for  the  mean  equation. 

Consider  two  fields  ^  ( sc , oj )  and  ip*  (x^ , &)  .  The  field  \p(x^fui) 
is  observed  at  point  _xt,  at  frequency  to,  while  the  field 
is  the  conjugate  of  the  field  observed  at  point  x^,  at  frequency 
to.  In  the  remainder  of  the  Appendix,  the  dependence  of  \ p  on 
xt  and  to,  and  of  on  x^  and  to  will  be  suppressed  for  conven¬ 
ience.  All  quantities  without  the  conjugate  sign  will  be  assumed 
to  be  functions  of  unhatted  variables,  while  conjugated  quantities 
will  be  functions. of  hatted  variables. 

The  total  field  t p  can  be  decomposed  into  two  parts:  the 
incident  field,  <p±,  and  the  various  scattered  fields  from  the 
individual  particles,  i.e. 

N  ... 

ip  =  \p.  +  £  ,  Lip.,  =  g  (A-l) 

j=l 

where  is  the  scattered  field  from  the  jt^1  particle  which 

is  related  to  the  transition  operator  Tj  of  the  jfc^  particle 
as. 

L*ij)  =  9^  -  Tj*(j)  (A-2 ) 

The.  quantity  in  (A-2)  is  the  field  at  the  location 

of  the  jth  particle  with  the  jth  particle  removed.  Substituting 
(A-2)  into  (A-l)  gives  the  following  equation 

N 

ip  =  +  Z  L^T.^*  ,  Lt{/.  *  g  (A-3) 

1  j  =  l  1 

Similarly,  the  equation  for  the  conjugate  field  ip*  ■  \p*(x£,0>) 
is  obtained  as 


(A- 4  ) 


I  *  ,  *  V"'  t"1*™*  J])  *  T*l*  _* 

P -  'Pi  +  E  L  ,  L  ip.  =  g 

1  j=l  3  1 

Now,  forming  the  product  P’P  *  and  averaging  yields 

N  ...  N  ... 

<\pip*>  =  ip.ip*  +  Ip.  E  L~  <T*ip{J,*>  +  E  L~  <T.ip[J,>ip* 
11  1  j=l  3  j=l  3  1 


“l*—*, (j)*  T* ,* 

t.i|>  J  =  g 


N  N 


E  E  L"1L"1*<T.T^(j)^(k)*> 

i=l  k=l  3  * 


(A-5) 


The  double  sum  appearing  in  the  above  equation  can  be  broken 
into  like  terms  (j=k)  and  unlike  terms  (jj^k)  as  follows 

<ipip*>  —  ip '  ip*  +  ip  _  V  L_1*<T*^(j)  *>  +  E  L~1<T.ip(j)>iP*  +  E 
11  1  j=l  3  j=l  3  1  j*l 

•  <T  .T*i{^3*  ip^3^  *>  +  E  E  L_1L“1*<T.T>(j)^(k)*>  (A-6  ) 

3  3  j=l  k=l  3  * 

To  obtain  an  equation  for  the  correlation  function  a  closure 
approximation  is  employed.  The  approximation  to  be  used  is 


i p'i'  z  <i p> 

^(j)^(k)*  ^  <(jJ><lp*>  f  j^k 


( A-7 ) 


This  is  a  generalization  of  Foldy's  closure  approximation 
used  to  derive  the  mean  equation.  Using  this  approximation 
in  (A-6)  gives 


IN  j  « 

<H*>  -  1> i<Pi  +4;  E  L”  *<tJ><i)*>  +  E  L"  <T ->«p>ip*. 

j  =  l  3  j=l  3  1 


•  •  1  •  1  h  it  a  • i<a^  ■  ^  •*  *  ^  1  #  it 

+  EL  L  <T  .T  .xipip*  >  +  E  El  <T,><ip>L  <T  .Xip  > 

j  =  1  3  3  j  =  l  k*l  ”  3  .  '  3 

j*k 


N  N 


.*-!*_*  , 


( A-8 ) 
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Since  the  particles  are  identically  distributed,  the  sta¬ 
tistics  for  all  particles  are  the  same,  i.e.. 


N 


£  <T • >  =  N<T> 

j=l  3 


N 


£  <T.T.>  =  N<TT*> 
j=l  J  D 


( A-9  ) 


{ A-10 ) 


Employing  these  simplifications  in  (A-8)  yields 
=  4,^4,*  +  NipiL_1*<T*><4>*>  +  NL~ 1  <T><ip>\p? 


+  NL~1L~1*<TT*><ip\p*>  +  {N2  -  N)  L_1  <T><4»>L~ 1  *<T*><ip*> 


( A-ll ) 


Notice  that  for  N»l,  N2  -  N  -  N*(l  -  J)  -  N2.  Using  this 
fact.  Equation  (A-ll)  can  be  put  into  the  following  form 

<<pip*>  zt  (ip^  +  NL~  *  <t><ip>)  •  +  NL-1*<T*><4/*>) 

+  NL~'L'”1*<TT*><4/4>*>  (A-12 ) 


Equation  (A-12)  can  be  simplified  substantially  by  making  use 
of  the  mean-wave  equation. 

The  mean-wave  equation  is  given  in  the  text  by  Equation 
(2-2-48)  as 

£  <ip>  =  g  ,  £  =  L  -  N<T>  (A-13 ) 

Multiply  (A-13)  by  from  the  left.  This  yields 

<i»  =  ! Jk  +  NL"’<T><i|»>  ( A-14  ) 

A  similar  result  for  the  conjugate  field  is  also  obtained 
by  inspection 

<tp*>  a  41*  +  NL~  *  *<T*x<J/*> 


( A-15 ) 


Substituting  Equations  (A-14)  and  (A-15)  into  Equation 
(A-12)  gives  the  desired  equation  for  the  space-frequency  cor¬ 
relation  function  <tptp*> 

<lplp*>  =  <  vp  >  <  *  >  +  NL”  1  L~  1  *<TT*><i|nj>*>  (A-16) 

Multiplying  Equation  (A-16)  by  ££*  from  the  left,  using 
the  fact  that  £  <ip>  =  g,  £*<i p*>  =  g*  and  ££*L  L  ~  1  yields 

££*<\p'p*>  -  N<TT*><t|)i p*>  -  gg*  (A-17) 

The  above  steps  can  be  justified  by  scaling  the  problems 
with  respect  to  the  characteristic  size  of  the  particle,  then 
N  is  replaced  by  the  fractional  volume  which  is  small. 

The  final  form  of  the  correlation  equation  as  it  appears 
in  Equation  (2-2-44)  is  obtained  by  writing  the  spatial  average 
of  TT*  out  explicitly  and  by  using  Equation  (2-2-38). 


APPENDIX  B 


Radon  Transform  of  the  Mean  Green* s  Function 

given  in 


(B-l) 

Employing  the  two-variables  perturbation  method,  the  mean  Green's 
function  G(x,x0)  is  expanded  as 

oo 

G (x  -  x0 ;x,x0 ; e)  =  £  G  n} (x  -  x0 ;x,x„) *en  (B-2J 

n=0 

where  the  slow  variables  are  defined  as 

X  =  ex  ,  x0  =  e  x o  ,  s  =  es  ,  p(s)  =  ep(s)  ,  e<<l  (B-3 ) 

Here  e  is  the  fractional  area  and  can  be  introduced  into  the 
Green's  function  equation  by  scaling  variables.  Setting  p  *  e p 
is  a  more  direct  but  less  physical  way  of  obtaining  the 
cprrect  e  dependence. 

The  Radon  transform  (  J  of  G^n^  (denoted  by  G^n^)  with 
respect  to  x  -  x0  is  defined  as 

RG(n)  s  G(n)  U,o;x,xJ 

=  y*dxG(n)  (x  -  x0;x,x0)  5  U  -  o*  (x  -  xa) )  (B-4  ) 

Thus  the  Radon  transform  of  the  mean  Green's  function  is 
found  to  be 

00 

RG  H  GU,o;x,x,)  =  £  G(n)  U,o;x,x0)  *en  (B-5) 

n=0 

In  this  Appendix  the  Radon  transform  of  the  mean  Green's 
function  is  approximated  to  the  first  order  only  as 


The  two-dimensional  mean  Green's  equation  is 
explicit  form  as  follows 

(V*  +k„)G{x,x0)  +  y*dsdx'p(s)t(x -s,x'  -s)G(x’,x0) 

=  -  6  (x  -  x„) 


After  inserting  Equation  ( B— 2 )  in>_o  (B-l)  and  equating 

o  ( °  ) 

coefficients  of  e  and  e'  ,  the  following  two  equations  in  G 
and  G ^  1  ^  are  obtained  respectively 

(V£  +  k„)G^°'  (x  -  Xo  ;x,xc  )  =-6(x-x0)  (B-7a) 

(V?  +  ko  )G  ^  1  ^  (x  -  x0  ;x,x0)  =*  h  (x  -  x„  ;x ,x0 )  (B-7b) 

where 

_  __  —  (  O  )  —  — 

h(x-x0;x,x0)  =  -2Vfc*VtG  (x-x0;x,x0) 

_  _  r  (o)  —  — 

-p(x)  /  dsdx't (x-s,x*-s)G  '  (x->:0;x,x0) 

V.  =  9  x  +  3  y  ;  V.  =  3  -  x  +  3-  y 
t  x—  y—  t  x  —  y  *- 

The  Radon  transform  of  Equation  (B-7a)  with  respect  to 
x  -  Xo  is  found  to  be 

(~-  +  ko)G(o)  (r,,o;x,Xo)  =  -6(C)  ,  -»<c<»  ( E  *10 ' 

(o)  _ _  (  o  )  _  _ 

Since  G  ( -  c,,  _o;XyX.° )  =  G  (  C»o;3C,x0 ) »  the  solution  of  (B-10) 

has;  the  following  form 

’  A(x;o)e^ot  , '  for  C  >  0. 

G ^ ° ^ ( C ,o;x ,x„ )  =  1  (B-ll) 

,  A(x;-o)e~jkoC  ,  for  C  <  0 

where  A(_x,o)  is  to  be  determined.  Notice  that  the  dependence 
of  A(>c;o)  upon  Xo  is  suppressed  for  convenience  and  the  solutions 
e  -*kof’  for  ^  >  0,  e^k°r'  for  5  <  0  have  been  excluded  since  there 
are  no  boundaries  in  free  space. 

/s  (  °  )  # 

The  continuity  of  G  at  r,  =  0  together  with  the  jump 

(  O  j 

condition  of  3  G  /3?;  at  c  =  0  imply  that 

A  ( x „  ;  -o )  =  -  ( B-12 ) 


A  ( x 0  ;  o)  * 


( B-8  ) 
( B-9 ) 


B-2 


The  Radon  transform  of  Equation  (B-7b)  with  respect  to 
x  -  x0  is  found  to  be 


+  k^G*  (C,o;x,x0)  =  h ( £ ,o;x,x0 ) 


(B-13 ) 


where  h(£,o;x,x.0)  is  the  Radon  transform  of  h  given  in  Equation 

A 

(B-8).  The  Radon  transform  h  is  computed  for  ;  >  0  as 

~  _  _  _  3ko? 

hU,o;x,x0)  =  A  ( x ; o )  e 


-  2iTjp(x)  J dst  (so,so)  f  (x;o,s)e-*S^  (B-14) 


where 


A(x;o)  =  -  2 jk0o*  VfcA(x;o) 


-  (2tt)  !p(x)  t (k0o,k0o) A(x;o) 


(B-15) 


In  Equations  (B-14)  and  (B-15),  t  is  the  Fourier  transform 
of  the  transition  kernal  t  of  the  scatterer,  and  the  function 
f(x;o,s)  is  defined  as 


f  ( x ;  o ,  s )  = 


A(x;o)  A (x; -o) 
s-k„  +  s  +  k„ 


(B-16 ) 


To  eliminate  secular  terms  from  the  pertubation  expansion 
of  Equation  (B-5),  the  following  secular  condition  is  required 


lim 

x  is  fixed 


G  ^ 1 ^  (C,o;x,x0) 


(B-17) 


This  secular  condition  implies  that 


f(x,o)  =  0 


(B-18 ) 


Inserting  Equations  (B-14),  (B-16)  and  (B-18)  into  (B-13.) 

and  defining  o  •  j<  =  jf  ,  o-x0  -  ifo,  the  following  differential 
equation  in  A (X;o)  for  z,  >  0  is  obtained 


B-3 


( B-19 ) 


dA ( c ; o)  2_2,  _ 

- — - '=  ~tz  pCx(^;o)}t  (k0o,k0o)A(C;o) 

d?  0 

Solving  Equation  (B-19)  with  condition  (B-12)  yields  the 
following 

_  i  jA$(C»Co»  o ) 

A(C;o)  =  -  23^-  e  (B-20) 

where 

A<J>(C,C0;o)  =  ~~  J  d? 'P  (xU' ;q)^t(k0o,k0q)  (B-21) 

Co 

Inserting  Equation  (B-20)  into  (B-ll)  and  by  Equation  (B-6), 
the  Radon  transform  of  the  mean  Green's  function  is  found  as 


RG  =  G(o)  U  ,o;x,x0)  = 


“  2jTc7  exp[jkoC  +  jA<HC,Co;o)]  ,  C>0 

— —  exp  [-jkQC  - jA(j>  (?,?„, --o)]  , 


( B-22 ) 


2jk0 


C  <  0 


APPENDIX  C 


Two-Dimensional  Transport  Equation 

Starting  from  the  correlation  equation,  the  two-dimensional 
transport  equation  can  be  derived  using  the  two- variable  pertur¬ 
bation  method  and  the  Radon  transform  technique. 


An  approximate  equation  of  the  Bethe-Salpeter  type  for  the 
correlation  function  , 

r**t'^t*  =  (^t)  >  (C-l) 

is  given  by 

££*T -J*dstp(st)T(st)T*(st)r  =  0  ( C-2 ) 

where  £  is  the  mean  field  operator  defined  in  Equation  (6-1-16) 
and  T(st-)  is  the  transition  operator  defined  in  Equation  (6-1-7). 
Equation  (C-2)  can  be  rewritten  as 

(£  -  £*)r -(£*”‘  -  £)  Jdstp(st)T(st)T*(st)r  =  0  (C-3) 

where 

(£  1'P)  (xfc)  =  j dx^  G  (xt,xp  ip  (x£)  ( C— 4  ) 

Here  G(sct,jc4)  is  the  mean  Green's  function.  It  satisfies 


£G(xt,x£)  =  Mxt  -  x£)  (C-5 )' 

Equation  (C-3)  is  now  written  in  a  more  explicit  form  as 

(V^  -  T (xt,xfc)+M(xt,xt)-M* (xt,xt)-N(xt,xt)+N* (xt,xfc)  =  0  (C-6) 

where 

M<xt,xt)  =  f  dxtldstp  (sfc)t(xt  -  St,xfcl  -  st)  T  (xu,xt)  ,(C-7) 

and 

N(xfc,xt)  = 


f^—tl 1 2 d^t 2 t p  (it°  G  (^t'5tl) 


(C-8 ) 


t^-tl  " -t'-t2  ” -t^  ^-tl  "  it'-t2  “  £t^  r  ^-t2'-t2^ 


For  small  fractional  volume  e ,  the  solution  of  Equation 

(C-6)  has  two  scales  of  variation  in  both  2it  ah<3  j<t*  A  problem 

of  this  type  can  be  treated  by  the  two- variable  perturbation 

, » 

technique.  The  small  parameter  e  is  introduced  into  Equation 
(C-6)  by  letting 

p(sfc)  =  ep(st)  ,  £t  =  est  ( C-9 ) 

where  it  has  been  assumed  that  p  is  a  function  of  the  slow  vari¬ 
able  £t  only.  This  assumption  implies  that  T,  G,  M  and  N  are 
functions  of  not  only  x^-,  x^  and  e  but  also  the  slow  variables 

A 

xj-  and  xt ;  thus 

A 

A(xt,xt;c)  =  A(xfc  '  (C-10) 

A  e  { r , G,M,N] 

The  function  A  is  now  expanded  in  a  power  series  in  E: 


/v  00 

A(*t'-^t;^t'£t;e*  =  (£t'£t?£t'^t^"  '  ( C— 11 ) 


A  e  { F ,G,M,N] 


Employing  the  chain  rule  the  and  operators  become 


where  Vt  and  Vt  are  the  del  operators  with  respect  to  the  slow 
variables  x^-  and  x^  respectively. 


equal  powers  of  e 
equations  for  the 


V2'  ■>  +  2eVfc  •  Vfc  +  e 2  v£ 


(C-12 ) 


V2  -*•  V2  +  2c V.  •  V.  +  e2V? 

z  z  z  z  z 


(C-13 ) 


Proceeding  v'ith  the  perturbation  analysis,  the  expansion 
given  in  Equatio^i  (C-ll)  is  used  in  Equation  (C-6).  Equating 

and  using  the  fact  that  =  n(°)  =  0  gives 

rf1).  The  first  two  of  these  are: 


(o) 


(V2  -  V2 )  r  (x.  ,x.  ;x.  ,x.  )  = 


0 


( C-14 ) 


m 


Q 


c 


c 


L* 


/  |  »  A,  ^  A  /  t  A 

(v£-v£)rl  (xt#xt;xt,xt)  +  2  (Vt*Vt-Vt*7t)r  K°}  (xt»xt;xt,xt) 

Ml  ^  (1)*  /y  _  (1)  ~  _ 

+  M  (xt,xt;xt,xt)  -  Mv  (*t'*t;-t'-t*  ” N 

/ 1  \  ^  _ 

+  N  t'— 1#— =  ® 

Equations  (C-14)  and  (C-15)  are  solved  using  the 
Radon  transformation  defined  as 


(C-15) 

double 


RA  *  A(^ ,6;o,o;xt,xt) 


dxt<3xt<5  (C  -  o-xt 


A 

)  6 (C  -  o-xt)A(xt,xt;xt,xt) 
A  e 


(C-16) 


Here  o  and  o  are  unit  vectors.  The  solution  of  Equation  (C-14) 
in  terms  of  the  Radon  transform  is  given  as 


?(o)  U>£;o,o;xt,xt)  =  IJ+  (xtfXt;£)eik  ^  ^ 


(C-17) 


where  J+  and  J~  are  still  to  be  determined.  The  Radon  Transform 
of  Equation  (C-15)  is 

(— - .  -2ik(o-7f +5-f)f  <") 

VH2  h2J  -t- 

-  fi<'>  -«<■>*  - S<‘>  ♦  S(,>*  (C-18) 

where  M^1),  1 )  are  the  Radon  transforms  of  and  N^1). 


To  insure  that  f  )  is  a  bounded  function  of  c  and  £  »  the 
right  hand  side  of  Equation  (C-18)  must  be  orthogonal  to  the 
homogeneous  solutions  of  Applying  this  so-called  secular 

condition  results  in  the  desired  transport  equation  for  the  two 
point  specific  intensity  J(xt»_xt»c>)  where 


C-3 


C 


J  (xt#xt;o)  ,  0  <  0  <  it 


J  (x.t,xt;o)  ,  ir  <  0  <_  2ir 


(019) 


Here  J+  and  J~  are  the  unknown  slowly  varying  coefficients 
appearing  in  Equation  (C-17).  The  transport  equation  for 
J(x^,^tJ£)  is: 

2iko*  (Vfc  +  Vfc)  J  fxt»xt;o)  +  (2it)  2p'{xt)  t  (ko,ko)  J(xt#xt;o)  -  (2ir)2p(xt) 


•t  (ko,ko)J(xt,xfc;o)  -  4i7r3p(xt)J  d9  •  1 1  (ko,ko/ )  |  2J (xfc ,xt ;o* ) 
+  _  ^  _  -  _ _  r2ir 


*A+(xt'*t;-)Ut~(-t  ~  3  “  4iif3p  (xt)f  d8 '  1 1  (ko,ko ' )  | 

J0 

—  —  +*  —  _  _ 

*J(xt#x  ;o')A  (x  ,x  ;o)U[o(x  -  x  ) ]  =  0 


(C-20) 


Here  t(Jc,k' )  is  the  Fourier  transform  of  t(x,x'  )  as  defined  in 
Equation  (6-1-12)  and 


+  —  — 

+  lA(^  (x.  , x.  jo) 

A  (*t»x;o)  =  e  ~z 


(C-21 ) 


with 


A^+  (xt,xt;o)  =  ± 


d?' P  (C  '  ;o) t (±ko;±ko)  (C-22) 


C  *  o-xfc 

The  unit  step  function  is  defined  by  U(C). 

1  ^ 

The  correlation  function  °  ( x ,  x ;x:,£)  is  found  by  taking 
the  inverse  Radon  transform  which  is  given  by 


rU)(: 


^  ^  TT  t  *  A  I  »  A 

~  4J  d2J  do(KKr(,,(C,C;o,o;xt,xt)  (C-23: 


0  0 


where 


C  =  o  •  x  ,  C  »  o  •  x 


(C-24) 


(C-25) 


(KK  T  '  '  )  (?,?  ;o,o;x,x)  = 


^  ^  {  fr  i-y 


(? fC»o»o» ’ 


with  pj*  being  the  principal  value  integral.  Substituting  Equation 
(C-17)  into  Equation  (C-25),  the  correlation  function  can  be 

/v 

obtained  in  terms  of  J(xt,xt;o)  as 


r(°)( 


iko* (x. -x. ) 


—  k2f  —  —  i  x ,  -x . 

x. ,x  ;x  ,x  )  =  -  I  do J  (x.  ,x.  ;o)e 

■  4  it  2  J  ~  “t  -t  - 


(C-26) 


In  summary,  the  correlation  function  is  approximated  by 
the  first  term  in  a  two  variable  expansion 


r(*t'^t*  =  r *°*(£t'-t;-t'-t> E°  + 


(C-27) 


where  r  (  °  J  (x.t'2Lt'2Lt  »2Lt )  is  9iven  in  terms  of  the  two-point  specific 
intensity,  J(xt,xt;o),  as  expressed  in  Equation  (C-26)  and  the 

A 

two  point  specific  intensity  J(xt,xt;o)  in  turn  satisfies  the 
generalized  transport  Equation  (C-20). 
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